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SUMMARY 


OB.TECTIVE 

The  main  objectives  of  this  research  program  are 
(1)  to  obtain  an ' analytical  solution  for  the  nonlinear 
response  and  ultimate  loads  of  two  hinged  reinforced 
concrete  circular  arches,  subjected  to  static  and  dynamic 
loads;  (2)  to  determine  experimentally  the  static  and 
dynamic  ultimate  loads  of  two-hinged  semicircular  reinforced 
concrete  arches  under  certain  typical  loading  conditions 
and  (3)  to  estimate  experimentally  as  well  as  analytically 
the  values  of  the  ductility  factor  )x  at  failure, 
being  defined  as  the  ratio  of  maximum  deflection  at  failure 
to  the  elastic  limit  deflection. 

The  nonlinearity  of  the  response  is  due  to  the 
nonlinear  stress-strain  curve  of  concrete  and  also  due  to 
the  effect  of  large  deflections  on  the  strains  and 
equilibrium  of  an  element  of  the  arch. 

SCOPE 

The  analytical  part  of  this  research  program 
includes  the  formulation  of  the  governing  equations  and 
boundary  conditions  for  both  static  and  dynamic  cases, 
taking  into  account  the  nonlinear  stress-strain  curve  of 
concrete  and  the  effects  of  large  deformations  on  the 
strains  and  equilibrium  of  the  arch  element.  The  governing 
equations  are  obviously  nonlinear  and  discontinuous  and 
are  solved  by  numerical  methods.  Programs  for  the  IBM  7090 
digital  computer  are  prepared  for  computing  the  response 
and  the  ultimate  loads,  using  these  numerical  methods. 

The  experimental  program  includes  the  testing 
of  small  scale  semi-circular  arches  to  determine  the 
static  and  dynamic  ultimate  loads  and  the  ductility  factor  }x 
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for  the  following  loading  conditions: 

1.  Uniformly  distributed  symmetric  radial  loads. 

2.  A  concentrated  load  at  the  crown, 

3.  Antisymmetric  concentrated  loads  at  quarter 
points . 

CONCLUSIONS 

The  experimental  and  analytical  studies  indicate 
that  the  approximate  conventional  theory  based  on  limit 
analysis  is  quite  adequate  to  predict  the  static  ultimate 
loads  of  the  underreinforced  arches.  The  dynamic  ultimate 
loads  for  compression  mode  loading  can  also  be  predicated 
by  approximate  theory  provided  that  an  appropriate  dynamic 
increase  factor  (based  on  the  increase  in  material  proper¬ 
ties)  is  used.  However,  it  is  not  clear  whether  a 
satisfactory  approximate  theory  can  be  developed  to  predict 
the  ultimate  dynamic  load  carrying  capacity  of  the  arch  , 
subjected  to  a  concentrated  dynamic  load  at  the  crown  or 
antisymmetric  concentrated  dynamic  loads  at  quarter  points. 

The  experimental  and  analytical  investigations 
also  indicate  that  the  natural  periods  of  the  arches  have 
a  significant  influence  on  the  dynamic  load  carrying 
capacity  of  the  arches  in  the  cases  of  a  concentrated 
load  at  the  crown  and  antisymmetric  concentrated  loads  .j 
at  quarter  points.  The  mode  of  failure  under  both  static 
and  dynamic  loads  is  quite  ductile  in  these  cases  as 
compared  to  the  compression  mode  loading. 
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CHAPTER  1 


INTRODaCTION 


1.1  OBJECTIVE 


e 


This  research  program  has  the  following  threefold 
objectives  s 

a)  To  obtain  an  analytical  solution  for  the  non¬ 
linear  response  and  the  ultimate  loads  of  two-hinged 
circular  reinforced  concrete  arches  under  static  and 
dynamic  loading.  The  nonlinearity  of  the  response  is 
obtained  by  including  the  effects  of  the  nonlinear  stress- 
strain  curve  of  concrete  under  compression  and  those  of 

# 

large  deflections  on  the  strains  and  the  equilibrium  of 
an  element  of  the  arch. 

b)  To  determine  experimentally  the  static  and  dynamic 
ultimate  loads  of  two-hinged  semicircular  reinforced  concret 
arches  under  certain  conditions  of  loading. 

c)  To  obtain  the  analytical  and  the  experimental  values 
of  the  ductility-factor  JJ^  at-  failure,  AX  being  defined 
as  the  ratio  of  the  maximum  deflection  at  failure  to  the 
elastic  limit  deflection. (Refer  to  §  2.6) 

1.2  PREVIOUS  WORK 

1.2.1  Analytical  Work; 

(1)* 

In  1932,  Cross  and  Morgan  summarized 

Superscripted  numbers  in  parentheses  refer  to 
references  given  In. the  bibliography. 


methods  of  analysis  and  design  of  reinforced  concrete  arches 
under  static  loading.  These  methods  were  based  on  the  con¬ 
ventional  linear  elastic  theory  for  reinforced  concrete 

structures.  In  19^1,  a  special  committee  of  the  American 

(2) 

Concrete  Institute  published  a  joint  report  on  plain  and 
reinforced  concrete  arches.  This  report  recommends  a 
numerical  procedure  to  account  for  the  moments  caused  by  the 
axial. thrust  and  the  deflection  of  the  arch.  It  also  suggests 
the,  use  of  Whitney's  stress  block  method  for  obtaining  the 
ultimate  strength  of  the  cross-section  of  the  arch.  In 
1953>  Onat  and  Prager  published  a  paper' concerning  limit 
analysis  of  arches  constructed  from  homogeneous  elastic- 
plastic  materials.  This  paper  proposes  a  theory  to  account 
for  the  reduction  of  the  plastic  moment  of  a  section  which 
is  subjected  to  both  a  moment  anS  a  thrust.  However,  since 
in  this  paper  the  properties  of  materials  are  considered  to 
be  Identical  in  compression  and  tension,  the  proposed  theory 
cannot  be  applied  to  reinforced  concrete  arches.  In  I960, 

Jain  published  a  paper^^^  on  ultimate  strength  of  reinforced 
concrete  arches.  In  this  paper  the  author  has  employed  a 
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bilinear  elastic-plastic  stress-strain  curve  £or  concrete 
and  obtained  the  ultimate  loads  making  use  of  an  iterative 

procedure  to  account  for  large  plastic  deformations  of  the 

# 

arch. 

( 

In  the  flfeld  of  dynamic  response.  Love  has 

presented  the  mathematical  formulation  of  the  vibration  of 

a  circular  elastic  ring  and  obtained  the  normal  functions 

of  free  vibrations.  In  I960,  Eppink  and  Veletsos  published  a 

paper on  dynamic  analysis  of  circular  elastic  arches.  Later 

(?) 

this  work  was  continued  to  include  inelastic  effects  The 

authors  have  considered  a  material  with  bilinear  stress- 
strain  curve  and  identical  properties  in  compression  and 
tension.  They  have  developed  a  set  of  governing  equations  to 
include  the  effects  of  large  deflections  and  have  used  a 
numerical  method  to  solve  these  equations. 

1.2.2  Experimental  Work: 

A  series  of  large  scale  static  tests  on  reinforced 

{ 19) 

concrete  arch  ribs  and  three  span  reinforced  concrete 

arch  bridges^ were  conducted  by  Wilson.  In  this  work 
the  measured  values  of  the  reactions  and  stresses  due  to 
unit  loads  were  found  to  be  in  close  agreement  with  the 
theoretical  values  based  on  linear  elastic  theory.  Jaln^^^ 
tested  two-hinged  reinforced  concrete  arches  to  failure  and 
showed  good  agreement  between  results  obtained  from  the 
tests  and  those  obtained  by  using  his  proposed  theory. 


In  the  field  of  dynamic  testing  Technical  Report 
No.  2-590  of  Corps  of  Engineers  discusses  the  results  of 
tests  conducted  on  underground  and  burled  fixed-end  and 
two-hinged  reinforced  concrete  arches.  This  test  program 
was  a  part  of  operation  PLUMBBOB  in  1957* 

1.3  SCOPE 

1.3*1  Analytical  Study; 

A  set  of  equations  and  boundary  conditions  governing 
the  behavior  of  a  reinforced  concrete  arch  are  formulated. 

The  equations  are  derived  from  the  conditions  of  equilibrium, 
strain-displacement  relations  and  force-strain  relations. 

The  equations  in  cases  of  the  static  and  the  dynamic  response 
differ  in  that  In  the  latter  case,  the  equilibrium  equations 
include  the  effects  of  inertial  forces  and  the  dynamic 
properties  of  both  concrete  and  steel  are  used  to  obtain  the 
force-strain  relations. 

A  simultaneous  solution  of  the  governing  equations 
yields  the  response  of  the  reinforced  concrete  arch.  The 
equations  being  nonlinear  and  discontinuous  they  are  solved 
by  numerical  methods.  Programs  for  the  IBM  7090  digital 
computer  are  prepared  for  computing  the  response  using  these 
methods.  The  ultimate  load  of  the  arch  is  obtained  as  a  load 


The  force-strain  relations  are  discontinuous  since  they  consist 
of  three  groups  (Pigs.  2.5b,  2.5c,  and  2, $6),  each  valid  for 
a  different  strain  distribution. 
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under  which  the  Tnaxiraum  compressive  strain  in  the  arch  reaches 


or  exceeds  the  ultimate  strain  e  ,  of  concrete.  The  small 

u 

reserve  strength  which  exists  in  the  arch  at  this  stage  is 
neglected.  Because  of  the  particular  definition  of  ultimate 
load  used  in  this  work,  the  analytical  methods  developed 
can  predict  the  ultimate  loads  of  only  the  following  types 
of  loadings:  m 

a)  Uniformly  distributed  compression  mode 

loading  #'  p. 

b)  Anti -symmetrically  distributed  deflection 
mode  loading 

and  c)  Symmetrical  and  antisymmetrical  concentrated 
loads  except  a  single  concentrated  load 
at  the  cro\TO.  This  aspect  is  discussed . in 
detail  in  §  2.1. 3.  -  ■ 

1.3.2  Experimental  Study: 

The  scope  of  the  experimental  program  included  the 
testing  of  small  scale  model  arches  and  the ' determination  of 
the  static  and  dynamic  ultimate  loads  and  the  ductility- 
factor  XX  at  failure  for  the  following  conditions  of  loading: 

a)  Uniformly  distributed  symmetric  radial  loads 

b)  A  concentrated  load  at  the  crown 
and  c)  Antisymmetric  concentrated  loads  at 

quarter  points. 

In  the  dynamic  cases,  a  failure  load  pulse  of  triangular  shape 
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having  a  rise  time  between  ten  and  twenty  milliseconds  was 
applied.  Dynamic  tests  were  also  conducted  under  partial 
loads  for  the  uniformly  distributed  symmetrical  loading  case. 

1.3.3  Comparison  of  Analytical  and  Experimental  Studies: 

Within  the  scope  outlined  above  a  comparison 
between  analytical  and  experimental  Investigations  is  obtained. 
Also,  the  s^tlc  and  dynamic  behaviors  of  two-hinged  reinforced 
concrete  arches,  are  compared. 


CHAPTER  2 


THEORETICAL  INVESTIGATION 

2.1  MATERIAL  PROPERTIES 

In  this  article  the  nonlinearity  of  the  stress-strain 
curve  of  concrete  under  compressive  loading  is  discussed. 

Also  the  effects  of  rapid  rates  of  straining  on  both  concrete 
and  reinforcing  steel  are  presented  herein. 

2.1.1  Concrete; 

a)  Static  Behavior:  A  typical  stress-strain 
curve  for  a  concrete  cylinder  under  static  compressive  loading 

/  O  \ 

is  shown  in  Fig.  2.1.  Hognestad'  '  suggested  that  Ritter's 

parabola  was  a  good  approximation  for  the  curve  up  to  ultimate 

stress  and  that  Inge  Lyee's.  equation  for  the  initial  tangent 

modulus  was  satisfactory,  provided  that  f!  in  that  equation 

c 

was  replaced  by  0.85  f •  Hognestad  assumed  the  descending 

c 

portion  of  the  curve  to  be  a  straight  line.  However,  the 
shape  and  the  extent  of  this  portion  of  the  curve  is  both 
uncertain  and  difficult  to  measure.  Hence,  Hognestad's 
expressions  are  used  to  describe  the  stress-strain  curve 
upto  failure.  The  expressions  ar6 

f"  =  0.85  C2-  >a') 

c  c 

// 

£  =  E  =  I800OOO  +  -<'^60^  t=si  .  C2.->b^ 

50C  C  C 
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"f 


C2lc) 


/ 

e 

c 


C 

qnci  -  0.0033 


C^L.Iol:) 


C2.I  e) 


where  f  -  Static  ultimate  stress  of 

®  concrete  obtained  from 

cylinder  tests 

e'  =  Strain  corresponding  to 

°  stress  f^ 


eu  =  Static  ultimate  strain  of 

concrete 

^soc  ~  Static  Initial  tangent 

modulus  of  concrete 


and  fc  and  ec 


Corresponding  stress  and 
strain  on  the  static  curve 


Strain  ^  QncJ  e 

C  dt 

Figure  2,1  Static  and  Dynamic  Stress-Strain 
curves  for  Concrete 
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The  tensile  stresses  in  concrete  under  static 
loading  are  neglected. 

b)  Dynamic  Behavior:  A  typical  stress-strain 

curve  for  a  concrete  cylinder  under  compressive  loading  with 

IQ) 

a  rapid  strain  rate  is  also  shown  in  Fig.  2.1.  Watstein 
compared  the  compressive  strengths  of  concrete  cylinders 
tested  at  rates  of  strain  varying  from  a  low  value  of 
10  ^  in./(  iia. )  (  sec  . )  to  a  very  high  value  of  10  in. /(  in. )  (  sec . 
His  work  indicates  that,  with  very  high  rates  of  strain,  the 
dynamic  ultimate  strength  can  be  much  greater  than  the  static 
ultimate  strength.  Yang,  et  al^^*^^  suggest  that  the  parameters 
associated  with  the  dynamic  stress-strain  curve  may  be 
obtained  by  using 

.  .  .  (2.2.  a')  . 


and 


E  .  =  £ 

doc  30C 


E 

c 


C2.2-c) 
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f 

where  =  Dynamic  ultimate  stress  of 

concrete  obtained  from  cylinder 
testa  conducted  at  rapid  strain 
rates 

I  ” 

e^^  =  Strain  corresponding  to  stress  f^^ 

E,  =  Dynamic  Initial  tangent  modulus 

°  of  concrete 

4^^  •  =  Dynamic  strain  rate  of  concrete 

e  =  Static  strain  rate  of  concrete 

c 

■  Making  use  of  relations  similar  to  those  used  for 
static  behavior  we  obtain  the  following  expressions ^  ‘  • 


f  =  0,8 S’!  .  o') 

dc  ' 


Stress-strain  curve  for  unloading  and  reloading  of 
concrete  Is  assumed  to  be  the  same  as  that  for  the  initial 
loading.  Also,  the  tensile  stresses  In  concrete  are  assumed 
to  be  negligible. 
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2.1.2  Reinforcing  Steel: 

a)  Static  Behavior?  The  stress-strain  curve  for 
steel  under  static  compressive  and  tensile  loading  is  shown 
in  Fig.  2.2.  The  following  expressions  are  used  to  describe 
the  stress-strain  relation. 


-m 

f  =  E  e 

s  s  & 

■for 

e  <r  e 

and 

11 

for 

e  >  e 

(2.4) 


=  Static  yield  stress  of  steel 

=  Static  yield  strain  of  steel 

=  Modulus  of  elasticity  of  steel 

=  Corresponding  stress  and  strain 
on  the  static  curve. 


Figure  2.2  Static  and  Dynamic  Stress-Strain  Curves 

for  Steel 
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b)  Dynamic  Behavior:  The  effect  of  rapid  strain 
rates  is  to  increase  the  yield  stress  of  steel  above  the 
values  obtained  under  static  loading.  The  modulus  of  elasticity 
remains  practically  unaffected.  Figure  2.^  of  Ref.^^^^  shows 
the  effect  of  strain  rate  on  the  dynamic  yield  stress  of 
steel.  It  can  be  seen  that  for  both  structural  and  inter¬ 
mediate  grades,  the  Increase  in  the  yield  stress  above  the 
static  value  varies  between  20  and  when  the  time  to  reach 
this  yield  stress  varies  between  0.001  and  0.1  seconds. 

Hence  in  this  work  the  dynamic  yield  stress  is  assumed  to  be 
30^  larger  than  the  static  yield  stress.  Behavior  of  steel 
in  tension  and  compression  is  assumed  to  be  identical. 

The  dynamic  stress-strain  curve  for  steel  is 
also  shown  in  Pig.  2.2.  The  following  expressions  describe 
the  stress-strain  relationship. 


The  stress-strain  curve  for  unloading  is 

assumed  to  be  linear  (Fig.  2.2),  the  slope  being  equal  to  E  . 

3 

The  reloading  curve  is  assumed  to  be  linear 


with  a  slope  equal  to  when  the  dynamic  strain  is  leas  than 
the  previous  maximum  dynamic  stralnj  however,  when  the  dynamic 
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strain  is  greater  than  the  previous  maximum  dynamic  strain, 
the  reloading  curve  is  assumed  to  be  indentical  with  the 
initial  loading  curve. 

2.1.3  Failure  Criterion: 

A  failure  criterion  based  on  excessive  compressive 

strain  in  concrete  has  been  used.  It  is  assumed  that  failure 

occurs  when  at  any  section  of  the  arch  the  combination  of 

the  thrust  and  the  moment  produces  a  compressive  strain 

which  exceeds  the  ultimate  strain  for  concrete,  e  .  Such 

u 

a  criterion  would  Indicate  failure  when  the  strain  in  the 
extreme  fibre  of  any  section  becomes  greater  than  e^. 

This  criterion  neglects  a  certain  reserve  strength 
in  the  structure  because  even  after  the  strain  in  the  extreme 
fibre  exceeds  e^,  the  inner  fibres  up  to  the  neutral  axis 
have  low  compressive  strains.  Further,  in  a  statically 
indeterminate  structure,  more  sections  than  one  have  to 
fall  before  the  structure  collapses.  If  the  load  distribu¬ 
tion  on  the  structure  is  such  that  the  failure  of  the  necessary 
number  of  sections  (the  number  depending  upon  the  degree  of 
indeterminancy)  occurs  at  a  load  which  differs  only  a  little 
from  the  load  at  which  the  first  section  fails,  the .failure 
criterion  used  here  would  be  adequate.  Such  lo.ad  distri¬ 
butions  on  a  two-hinged  circular  arch  are, 

a)  uniformly  distributed  symmetrical  loads 

b)  uniformly  distributed  anti symmetrical  loads 
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and  c)  symmetrical  and  anti symmetrical  concentrated 
loads,  except  a  concentrated  load  at  the  crown. 

This  failure  criterion  Is  inadequate  when  a  single 
concentrated  load  la  considered  because  when  the  maximum 
compressive  strain  at  the  section  under  the  load  exceeds  e^, 
the  other  sections  of  the  arch  have  low  compressive  strains 
and  hence  the  arch  has  a  considerable  reserve  strength. 

The  failure  criterion  thus  gives  a  lower  bound  of 
the  ultimate  load  for  the  above-mentioned  distributions. 
However,  for  design  purposes  this. may  be  a  realistic  limit. 
2.2  GEOMETRY  OP  THE  ARCH 


The  two -hinged  arch  under  consideration  is  qhown  in 


Fig. 2. 3*  The  geometry  of  the  arch  is  described  by  the 
following! 


r,  =  Mean  radius  of  the  arch 
JZS  =  Half  the  central  angle 
b  =  Width  of  the  cross  section 


t  =  Depth  of  the  cross  section 

and  0  =  Angle  subtended  at  the  center  by  an  arc 

between  the  left  support  and  any  point  on 
the  arch. 


Plgupe  2.3  Geometry  of  the  Arch 
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2.3  STATIC  RESPONSE 


2.3.1  Governing  Equations: 

a)  Equations  of  Equilibrium:  Pig.  2.1).  shows  the 


geometry  of  deformation  of  an  element  of  the  arch  and  the 
forces  acting  on  it. 


Forces  and  Displacements  are  Shown  In 
Positive  Sensre 

Figure  2.4  Forces  Acting  on  a  Deformed 
Element 


Considering  tiie  equilibrium  of  th.e  forces  in 
r  and  9  directions  and  the  moment  equilibrium,  the  following 
equations  are  obtained; 

sm  (^)  _  2.Q  5\».Q4|L)ca5CAci0) 

di©  cosrdSV  cos  CAd0)  _  dd  d©  cosC^)  s'tnCAci©)  =  O 

o(0  3.  v'  a© 

C2.6) 

b  r  d©  +  stn  (^Ad©)  _  2.Q  Cos  CAcl©) 

—  dQ  cos(^d|^')  co3  CAd0)  —  d  0  C.05  s''v>CAd0^  — O 

d©  ^  d0  ^  1  ^ 

C2.V) 

«"=>  ^ -<5-; f  ^  =o  C2.S) 

where  p  =  Radial  component  of  load  per  unit 

^  length  of  the  arc 


P9 

— 

Tangential  component  of  load  per 
unit  length  of  the  arc 

d9 

= 

Angle  subtended  at  the  center  by 
the  element  before  deformation 

Ad0 

= 

Change  in  d0  due  to  deformation 

N 

= 

Axial  force  on  the  section 

M 

= 

Bending  moment  on  the  section 

Q 

= 

Shear  force  on  the  section 
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b)  Strain-displacement  Relations:  Assuming  that 
the  normals  to  the  middle  surface  remain  normal  after  the 
deformation  and  that  the  shear  strains  are  negligible, 
expressions  for  the  strain  in  the  middle  surface,  e^  and  the 
change  in  curvature  %  are  obtained  as  follows  (Pig.  2.I4.): 

The  length  and  the  curvature  of  the  undeformed 

element  are 


<is  =  h  cl©  •  dQ  _  JL 

®  ~  V- 


Ca) 


The  squsire  of  the  length  of  the  deformed 


element  is 


(ds  +Aols')  =  (  —  a©  -Y  CdliJ  +  ~  Cb) 


where  u  “  Tangential  displacement  of  the 

element,  measured  positive  when,  the 
element  moves  clockwise 

and  w  =  Radial  displacement  of  the  element, 

measured  positive  when  the  element 
moves  toward  the  center. 

Since  the  compressive  strain  in  the  middle 


surface  is 


we  get 


C'-e 


-♦-Ads') 


CO 
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From  equations  (b)  and  (c),  neglecting  the 


Y* 

I 


1  ,  du. 

do 


Hence,  the  change  in  curvature. 


a. 


C2.\o) 


c)  Force-strain  Relations:  Pig.  2., 5a  shows  a 
symmetrically  reinforced  rectangular  concrete  section  of  an 
arch,  acted  upon  by  a  positive  bending  moment  M  and  a  positive 
thrust  N.  Depending  upon  the  relative  magnitudes  of  M  and 
N,  the  strain  distribution  across  the  section  will  be  as 
shown  in  either  Pig.  2.5b  or  Pig.  2.5c.  It  is  assumed  that 
the  ratio  of  r^/t  is  large  enough  to  neglect  the  nonlinearity 
of  the  strain  distribution  across  the  depth  of  section.  If 
the  thrust  N  is  negative,  it  is  possible  that  the  entire 
section  will  be  in  tension.  (Fig.  2.5d) 


Cast  1 

Ca*e  2. 

C 

(a) 

CV.) 

Cc) 

(A) 

Figure 

2.5  Arch  Section  and 

Possible 

Strain  Distributions 
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Noting  that  positive  bending  tnotnent  causes 
a  decrease  in  the  initial  curvature  of  the  arch,  we  have 


e 

o 


2 


ancA  ^4) 

t 


/  C2.U) 


J 


Using  the  stress-strain  properties  of  concrete 
and  steel  described  in  §  2.1.1a  and  2.1.2a  respectively,  the 
following  relations  between  M,  N,  e^,  e^,  and  e^  are  obtained. 

Case  1  -  Section  completely  in  compression 

i .  e  . ,  e  <L  €  and  £  ^  O 
4  ^  «  1  ^ 


fbt  if 


E  ,  A 

-^—7  -be)  -7  — ^ 

//  ^  2,  J  Ce-e') 


C2.\2) 


^  » 


p'bt*  4r  '  t 

c  c  4  ' 


Ce.ej 


2.  Ce  -  e  > 

4  » 


.  A. 
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where 


=  Area  of  tension  steel 


=  Area  of  compression  steel 

/A  +  a\ 

(  s  s) 


=  A 


bt 


n 


Modular  ratio,  E  /E 

'  s'  sc 


sc 


=  Secant  modulus  of  concrete 


S  e' 


/  i 


6  = 


Z  e' 


A  ^ 

Ce  -  e  ) 

^  A _ I  ' 


A  e 


/  X 


Case  2  -  Section  partly  in  compression  and  partly 

in  tension. 

i  •  e  • »  e  <  e  and  e  <  O 

■4  ^  U  I 


N 

C 


E.  h 


A. 


2f  - 

c 


ce  -e) 

^  I 


(2,14) 


and 


_  ,  5L  Te  («-A )  -  e  "1  +  —  i 

-  t  L  ^  d  Ce-e^ 


-4  V  A 
2Ce  -€)^  ^ 

4  V 


C2..\5) 
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The  steel  strain  e.  or  e,  in  the  equations 

«  3 

(2.12)  through  (2.15)  is  replaced  by  e  when  either  the  bottom 

y 

or  the  top  steel  yields.  . 


’  Case  3  -  Section  completely  in  tension 
l*e.  e  40  and  6  40 


o 


C2-\e) 


C2.I7) 


It  is  evident  that  in  order  to  satisfy  both 
(2.16)  and  (2.17)  simultaneously,  the  tensile  strains  -e^ 
and  -e,  must  be  below  the  tensile  vield  strain  ^e^. 

d)  Boundary  Conditions:  The  boundary  conditions 
which  govern  the  solution  of  the  two-hinged  arch  problem  are 

At  0  =  0  and  Q  =  2^ 

O 

u  =  w  =  M  =  0  (2.18) 
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2.3*2  Nondimensionalization  of  the  Equations; 

The  governing  equations  (2.6)  through  (2.18) 


are  converted  into  nondimens ional  form  using  the  following  no¬ 
tation: 


P  dl©  Cosec  ^  |^s>nCAd£>)N  -  CosCA^e^Q 

O 

-  k  die  doi  (^)  (“cos  CAd6)-^  +  si^CAae)  ^  1 

^  ^  ^  L  ao  de  J 


=  o 


C2.2.0 


d© 


A 


cos 


(x)^ 


i- 


c:ie 

3.t 


^ -  /■  <^0  \  aw 

S»r>  f  -^  )  — :r- 

^  ^  d© 


=  O 


C  2  .  X  a") 


Strain-displacement  relations: 


€  =  W  - 


au  „  1  ^  _ IL  i£W 

e  )  <p^ 


do  ^  \9Sd0 


Ca-xs) 


d  W 

d©"- 


-L  ^ 

^  d© 

o 


(2.24) 


Force-strain  relations; 

^  Ce  -te  ^ 

€  —  I  J 

°  % 


— .  y» 

Qnd  %  =_a.re-e') 


(2.26-) 
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Case 


N  = 


Case 


N  = 


M  = 


1  : 


If 


6:^6  and 

^  u 


e  >o 


A. 


LA.  .4.  (v-1;)  (e  -e  )  -h  A 


Ce  -e  ) 

A  I 


re  +  e  ') 

''I  4 '' 


a  Ce  -e;) 

A  > 


f,  K 


1  Ce  -e  f  ^ 

A  < 


C2.2-0) 


C2.27) 


C  2. 2S) 


C2.2<^) 
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Case  3  : 


If  e  O  and  6^0 

H  =  - ^  -l-e  )  CX.3o) 

z?"  ^  ^  i 

C 

t  ^  ^ 

c, 

and  boundary  conditions: 

At  0  =  0  and  0  =  2 

U  =  W  =  M  =  0  (2.32) 

2.3*3  Finite-difference  Formulation; 

In  order  to  obtain  a  numerical  solution,  the 
differential  equations  (2.20)  through  (2. 24 )  are  converted 
into  difference  equations.  The  form  of  the  difference 
equations  chosen  is  suitable  for  the  numerical  method  described 
in  the  next  article.  The  difference  equations  represent 
a  discrete  system  consisting  of  (m-l)  segments  denoted  as 
etc.  The  ends  of  the  segments  are  denoted 
as  i-1,  i,  i+1,  etc.  (Fig.  2.6).  Each  of  the  segments 
subtends  a  nondiraensional  angle  A0  at  the  center  of  the  arch. 
The  unknowns  M,  R,  Q,,  e^,  qt  ,  U  and  W  at  each  section 
are  denoted  by  subscripting  them,  N^,  etc.  The  loads 

and  Pg  on  each  segment  are  denoted  as  ^j,y  ^Qy 
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Figure  2.6  The  Discrete  System 
Prom  equation  (d)  of  §  2.3.1b, 
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Using  forward  differences, 
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Using  the  averages  of  the  values  of  N 
and  Q  at  Sections  i  and  i+1  and  the  forward  differences 
for  the  derivatives  of  N,  Q  and  M,  the  differential 
equations  (2.20),  (2.21)  and  (2.22)  are  transformed 
into  the  following  difference  equations. 
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The  strain-displacement  relations  (2.23) 


and  (2.24)  are  used  to  obtain  difference  equations  which 
relate  displacements  at  sections  i  and  i+1. 
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The  finite-difference  approach  used  above 
was  chosen  as  the  most  suitable  to  obtain  a  good  numerical 
solution  after  making  a  preliminary  study  of  various  possible 
approaches . 

2.3.U  Numerical  Method  of  Solution: 

A  two-hinged  circular  arch  has  a  degree  of 
static  indeterminancy  equal  to  one.  The  structure  is 
analysed  by  treating  the  reaction  at  the  right  hand  support 
as  the  redundant  (Pig.  2.7).  The  magnitude  of  the 
redundant  is  determined  by  first  assuming  a  certain 
value  and  then  refining  it  by  successive  iterations 
until  all  the  governing  equations  are  adequately  sat¬ 
isfied. 


Figure  2.7  Arch  with  External  Reactions 

The  method  can  be  described  as  follows: 

Trial  1  - 

(a)  The  arch  is  divided  in  (m-1)  segments 
and  the  nondlmensionallzed  radial  and  tangential  components 
of  the  load,  distributed  on  each  segment,  are  obtained. 
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(b)  An  initial  nondimensional  value  is 
assmed  for  the  horizontal  reaction  at  the  right  support. 

A  good  initial  value  can  be  obtained  from  the  linear 
elastic  analysis  of  the  arch. 

(c)  Using  equations  of  static  equilibrium  for  the 

entire  structure,  the  nondimensionalized  values  of  the 

other  three  reactions,  viz.  H  ,  V  ,  and  7  are  obtained. 

A  A  B 

(d)  Equations  (2.37)  are  now  evaluated  assuming 
(ActO),  =  0.  Making  use  of  V  ,  H  ,  M  and  the  constants 

J  A  A  A 


c,  ,  c  ,  ...,  equations  (2.34),  (2.35)  and  (2.36)  are 

A  3 

solved  for  (m-2)  segments.  Thus  the  forces  N  and  Q  and 
the  moments  M  are  known  at  all  the  m  sections. 

(e)  At  all  the  m  sections,  the  values  of  concrete 
strains  at  the  top  and  the  bottom  of  the  section  are 
calculated  assuming  that  the  concrete  stress-strain  curve 
is  linear  in  compression  and  that  the  tensile  stresses  in 
concrete  are  negligible.  The  considerations  used  for  this 
purpose  are  given  in  Appendix  I. 

(f)  Using  the  strains  obtained  above,  each  of  the 
m  sections  is  classified  as  either  dase  1  (Pig.  2.5b), 

Case  2  (Fig.  2.5c)  or  Case  3  (Fig.  2.5d).  Governing 
equations  (2.26)  and  (2.27),  (2.28)  and  (2.29)  or  (2.30) 
and  (2.31)  are  chosen  for  each  section  in  accordance  with 
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the  above  classification.  The  first  two  pairs  of  the 
equations  can  be  solved  by  the  »rrewton-Raphson  iteration’ 
method'  .  ( If  M  is  negative  the  terms  e^and  e^  are 

interchanged  and  the  same  equations  are  employed  along 
with  the  absolute  value  of  M. )  This  method  as  applied  to 
the  present  problem  is  explained  in  Appendix  II.  The 
strain  values  obtained  in  step  (e)  are  used  as  initial  values 
to  start  the  iterations.  The  last  pair  of  equations  being 
linear,  the  strains  obtained  in  step  (e)  satisfy  these 
equations.  Thus  the  solution  of  the  proper  pair  of 
equations  gives  strains  e^  and  e^^  for  all  the  m  sections, 

(g)  Using  the  concrete  strains  e^  and  e^ ,  the 
steel  strains  e^  and  e^are  obtained.  If  |  e^ji  or  le^i 

or  both  are  greater  than  |e  |  ,  the  governing  force-strain 

<J 

relations  are  altered  as  explained  in  2.3.1c  and  the 
nonlinear  equations  are  once  again  solved  by  the  ’Newton- 
Raphson  iteration’  method.  The  initial  values  for  starting 
the  iterations  are  the  solutions  obtained  in  step  (,f ) , 

The  values  of  e  and  e  thus  obtained  take  into  account  the 

4  1 

yielding  of  steel. 

(h)  The  quantities  and  ot  are  obtained  using 
equations  (2.2$)  for  all  the  m  sections. 

(i)  Using  equations  (2.38a)  through  (2.3Bf), 

>  Uj  ,  w/'(i  =  l,...,m)  and  (1  =  2,...,m)  are 

I 

evaluated  in  terms  of  one  unknown  viz.  W .  Since,  in  the 

f  * 

first  trial,  values  of  Ui  h  and  W^'s  are  not  known. 
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eoj^^s  in  equation  (2*38a)  are  assumed  to  be  equal  to  e^^'s. 

This  amounts  to  neglecting  the  n<.'  linear  terms  In  the  strain- 

displacement  relations  for  this  cycle  of  trial  1.  Eauatlons 

(2.39)  and  (2.40)  are  now  written  for  (m-1)  sections  to 

obtain  a  set  of  (2m-2)  equations  involving  (2m+l)  unknowns 

viz.  m  Uj^'s,  W^'s  and  W^.  From  equation  (2.32) 

we  know  that  U  =  W  =  U  =0.  Thus  the  unknowns  are  reduced 
1  1  m 

to  (2m-2).  On  solving  these  equations  simultaneously  the 
values  of  deflections  U  and  W  are  obtained  for  all  the 
m  sections. 

(j)  Making  use  of  equations  (2.33)  £uad  (2.38a) 
(^'^®)j,  (j  =  !»•••»  m-2)  and  e^j^  (i  =  1,...,  m)  are  calcu¬ 
lated  and  steps  (d)  through  (1)  are  repeated  using  the 
newly  calcualted  values  of  (  Ad©  ) and  ®oj*  The  deflec¬ 
tions  U' 3  and  W*3  thus  calculated  Include  large  deflection 
effects.  These  deflections  are  found  to  be  satisfactory  and 
a  second  repetition  of  this  step  is  not  needed. 

Trial  2  - 

(a)  If  the  value  of  Wj„  obtained  In  Trial  1  -  (  j), 
does  not  satisfy  the  condition  »  0  (equation  2.32)*, 


*In  the  present  method  a  tolerance  of  (  -  I/500  Inch  <^win  ^ 

1/500  Inch)  is  considered  to  be  satisfactory. 
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the  assumed  value  is  changed  by  a  small  percentage 

to  a  new  value  • 

(b)  Steps  (c)  through  (j)  of  Trial  1  are  now 

repeated,  the  only  difference  being  that  the  values  of 
(Aol9)j  in  step  (d)  [^equation  (2.33)3  and  those  of  e^^ 
in  step  (i)  [equation  (2.38a)3  are  obtained  by  using  the 
values  of  and  ' s  calculated  in  step  (j)  of  trial  1. 

(c)  If  the’  most  recent  value  of  the 

value  of  W  is  modified  once  again.  The  new  value  W 

QX  D  6 

is  obtained  by  extrapolating  linearly  on  the  basis  of  , 

W  and  the  deflections  ’s  associated  with  each 
61  m 

respectively. 

Using  H-- ,  steps  (c)  through  (j)  of 

P  3 

Trial  1  are  repeated  and  the  deflection  v/^  is  checked. 

This  procedure  is  continued  until  the  deflection 
falls  within  the  tolerance  limits, 

2.3.5  Calculation  of  Ultimate  Load: 

The  ultimate  load  carrying  capacity  of  an  arch 
for  a  particular  distribution  of  loading  can  be  obtained 
by  using  the  method  described  in  the  previous  article.  The 
procedure  consists  of  analysing  the  arch  for  Increasing 
values  of  load  until  the  ultimate  is  obtained.  The  ultimate 

load  is  assumed  as  that  which  causes  the  maximum  compressive 
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strain  in  the  arch  to  exceed  e^,  the  ultimate  concrete 
strain,  (  §  2.1.3) . 

2.3.6  Digital  Computer  Program: 

A  computer  program  prepared  to  handle  the  extensive 
calculations  involved  in  obtaining  a  numerical  solution  by 
the  above  method  is  presented  in  Appendix  III.  The  program 
consists  of  a  main  routine  and  six  subroutines.  The  flow¬ 
chart  of  the  program  is  also  presented  in  Appendix  III. 

2.3.7  Convergence: 

a)  Number  of  Discrete  Segments: 

Semicirciilar  arches  were  divided  into  16,  20, 

2i|  and  28  segments  and  analysed  for  various  load  distri¬ 
butions.  It  was  found  that  the  results  obtained  from  the 
last  three  cases  were  essentially  similar  and  a  choice  of 
20  segments  (each  subtending  an  angle  of  9°  at  the  center 
of  the  arch)  was  made  to  approximate  a  semicircular  arch, 

b)  Convergence  of  Iterations: 

In  the  numerical  method  described  in§2/3.-4> 
two  different  Iterative  procedures  are  used.  The  first  one 
concerns  the  solutions  by  means  of  'Newton-Raphson  iteration* 
[steps  (f)  and  (g)  of  Trial  l3.  Certain  difficulties  ex¬ 
perienced  in  obtaining  convergence  with  this'  iteration  have 
been  explained  in  Appendix  II  (  SA2.3). 

The  second  iteration  involves  successive  choice 
of  the  horizontal  reaction  until  sufficient  convergence  is 
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obtained  in  as  few  as  fo'ur  to  six  cycles  of  iteration. 
However,  for  certain  load  distributions,  as  the  load 
.approaches  the  ultimate  small  changes  in  the  assumed  value 
of  the  redundant  cause  large  changes  in  displacement  and 
it  becomes  difficult  to  obtain  the  convergence.  Consider, 
for  example,  a  semicircular  arch  under  the  action  of  a 
uniformly  distributed  antisymmetrical  loading  (Fig.  2.8a). 


Figure  2.8  Semicircular  Arch  under  a  Unlforln 

Antlayimnetrlcal  Load 

For  points  below  'i.  on  the  load-deflection 
curve  in  Pig.  2.8b  the  convergence  is  satisfactory.  When 
the  load  is  Increased  to  that  at  point  X,  the  combination 
of  moment  M  and  thrust  N  for  a  certain  trial  value  of  the 


horizontal  reaction  say  causes  the  maximum  compressive 

strain  at  sections  K  and  L  (quarter  points)  to  exceed  e^. 
That  is,  the  strains  in  many  fibres  at  these  sections 
correspond  to  the  drooping  part  of  the  concrete  stress- 
strain  curve.  V/ith  such  large  strains,  the  change  in 
curvature  at  the  sections  K  and  L  is  large  and  its  con¬ 
tribution  to  the  end  deflection  W  is  large.  As  the  value 
of  the  horizontal  reaction  is  changed,  the  moments  and  the 
thrusts  on  the  sections  change.  These  changes,  though 
small,  may  be  sufficient  to  cause  a  change  in  the  strain 
distribution  at  either  K  or  L  or  both,  such  that  the 
maximum  compressive  strain  is  less  that  e^.  Such  a  drastic 
change  in  strains  can  occur  because  the  concrete  stress- 
strain  curve  has  been  assumed  to  be  parabolic  with  a 
drooping  part,  the  slope  of  which  increases  very  rapidly. 
The  net  effect  of  this  large  change  in  the  strains  at 
K  or  L  is  to  affect  the  curvature  at  these  points  con¬ 
siderably  and  consequently  to  change  the  end  deflect!  on 
Wjjj  by  a  large  amount.  On  account  of  such  a  sensitivity  of 
the  deflection  to  the  changes  in  the  value  of  the 
redundant,  the  solution  tends  to  oscillate  or  even 
diverge.  However,  it  is  seen  from  Fig.  2.8b  that  such 
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a  difficulty  in.  convergence  occurs  only  in  the  vicinity  of 
X,  where  the  load-deflection  curve  is.  almost  horizontal. 
Therefore,  it  is  reasonable  to  take  the  load  at  X  as 
the  ultimate  load  of  the  arch,  and  a  good  convergence 
is  not  necessary.  The  method  thus  yields  a  close 
approximation  of  the  value  of  the  ultimate  load,  but 
is  unable  to  predict  the  load-deflection  curve  beyond  X. 
2.4  DYNAMIC  RESPONSE 

2.4.1  Governing  Equations: 

The  forces  actireg  on  an  element  of  the  arch 
deformed  under  the  action  of  time-dependant  loads  are 
shown  in  Pig.  2.9.  In  addition  to  the  external  loads 
and  the  internal  forces,  radial  and  tangential  inertia 
forces  are  shown  to  act  on  the  element.  The  rotational 
inertia  has  been  neglected.  The  geometry  of  deforma¬ 
tion  of  the  arch  element  is  similar  to  that  in  the 
static  case  (Pig.  2.4).  Also,  the  force-strain 
relation  now  depend  on  the  dynamic  stress-strain 
properties  of  both  concrete  and  steel. 

a)  Equations  of  Equilibrium:  By  con¬ 
sidering  the  equilibrium  of  forces  in  r  and  0  directions 
and  the  moment  equilibrium,  the  following  equations  of 
dynamic  equilibrium  of  the  arch  element  are  obtained; 
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Figure  2.9  Dynamic  Forces  Acting 

on  a 

Deformed  Element 
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v/here,  in  addition  to  the  notations  used  in  §  2,3.1,  the 
density  of  reinforced  concrete  is  denoted  by  P  and  the  time 
by  t. 


b)  Strain-displacement  Relations;  Since 
the  geometry  of  the  deformed  element  is  similar  to  the 
one  in  the  static  case,  the  strain-displacement  relations 
are  similar  to  equations  (2.9)  and (2.10)  except  that  the 
total  derivatives  are  replaced  by  partial  derivatives; 
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Also  equation  (d)  of  52.3. Ih  becomes 


Ade  =  \  ■  'bn  de 
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c)  Force-strain  Relations;  The  force- 
strain  relations  also  are  similar  to  those  in  the  static 
case  except  that  the  dynamic  properties  of  both  concrete 
and  steel  are  to  be  used.  That  is,  in  equations  (2,11) 
through  (2.1?)  f|(  and  e^  are  to  be  replaced  by  f^^  and 
e^^  respectively.  The  dynamic  force-strain  relations  are 
given  in  a  nondimensional  form  in  §  2.4.2. 

d)  Boundary  Conditions*  The  boundary 
conditions  governing  the  dynamic  case  are  the  same  as  in 
the  static  case,  viz. 

At  0  =  0  and  0  =  20 

'^O 

u  =  w  =  M  =  0  (2.47a) 

e)  Initial  Conditions;.  The  initial 
conditions  specify  that  the  arch  be  undeformed  at  the 
start  of  the  response,  l.e.  at  time  t  =  0, 
u=w=M=N=Q=0  everywhere  on  the  arch  (2.47b) 

2.4.2  Nondimensionalization  of  the  Equations; 

For  the  purpose  of  converting  the  governing 
equations  into  nondimensional  form,  the  following  notation 
is  used  in  addition  to  that  defined  in  §2.3.2  [equations 
(2.19  )D, 
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where  T  =  NatiJral  period  of  the  arch. 

Hence,  the  nondlraensional  governing  equations  are. 
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Strain-displacement  relations: 
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Force-strain' relations: 

As  stated  earlier,  the  force-strain  relations  are 

obtained  by  the  substitution  of  f^^  and  e^^  in  place  of 

f”  and  e'  respectively  in  equations  (2,11)  through  (2.1?). 
c  c 

The  form  of  equation  (2,11)  is  altered  to  suit  the  numerical 
method  described  in  2. 4. 4*  Thus,  the  force-strain  relations 
are. 
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Case  2- 
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=  — ^  --  (e  -e  ) 

Af  ^  ^  ^ 


C2--S1) 


■k$- 


In.  dases  1  and  2,  the  -steel  strain  e^  or 


e^  is  replaced  by  e^^  when  either  the  bottom  or  the 
top  steel  yields.  However,  in  Case  3,  both  the  tensile 
strains  -eg  and  -Sg  need  to  be  below  the  tensile 
yield  strain  order  to  satisfy  equations  (2.60) 

and  (2.61)  simultaneously. 

The  boundary  conditions  in  nondimenslonal 

form  are. 


At 


0=0  and 


0=2 


U  =  W  =  M  =  O' 


(2.62a) 


The  Initial  conditions  are, 

at  t  =  0,  tr=W=M=ij  =  Q=  0  everywhere  on  the  arch 

(2.62b) 

2,li,3  Finite-difference  Formulation: 

In  order  to  solve  the  governing  equations 
Using  numerical  techniques,  the  differential  equations 
(2.2|9)  through  (2.54)  ®re  converted  into  difference 
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equations,  suitable  for  the  method  presented  in  2.4.4* 

The  formulation  of  the  difference  equations  is  accomplished 
by  using  a  discrete  system  shown  in  Pig.  2.6,  The  notation 
followed  in §2.3.3  denoting  the  discrete  sections  by 
i-1,  i,  ifl  etc.,  and  the  discrete  segments  by  j-1,  3, 

3+1  etc.,  is  also  used  in  the  present  case.  In  addition, 
successive  time  intervals  are  denoted  by  k-1,  k,  k+1  etc. 
Tlierefore,  each  unknown  carries  two  subscripts  - 
either  subscript  i  or  j  depending  upon  whether  the 
unknown  pertains  to  a  section  or  to  a  segment  and  the 
subscript  k  specifying  the  time  interval. 

The  difference  equations  obtained  are  as 

follows ; 

Using  forward  differences,  equation  (2.54) 

becomes , 


(Adey 


(w  _zw.  +  w.  3 

i.+a,k  L  ,v. 

f  A© 
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By  using  central  differences  for  sections 
away  from  the  boundaries  and  forward  differences  for 
boundary  sections,  equations  (2,^2)  and  (2,53)  become. 


^  ^  =  w. 


_  -  Ui-,,0 

T-f  A© 


6  (.</>  A0)^  i  A© 


C2.C4) 


ci  =  - 


e  -  _ 

°  ~7  _ 

Ae 


2  A§)^ 


£  —  ^rn-ijk  ,  "Wm-i.k 

o  Tr\  y.  IT 

"  0.  C<jl>  A©) 


X.  =  e  .  i- 

Cjk  o  i^w, 


CjiS  h^f 


8C5z5  a©)"^ 

o 


t  k  oik  -L 

o-Cp 

o 

% 

_  p  _  ^^''^-1  JV. 

*n  k  "  k  _  1 

'  '  2  Cjl5  A©) 

o 


(2-9^  A©)  \  C2-^0 


(  L  =  2.^--- 


Using  the  notation  defined  by  equation  (2.37), 
equations  (2.49),  (2.50)  and  (2.5l)  are  converted  into 
difference  equations.  The  quantities  N  and  Q  are  sub¬ 
stituted  by  their  average  values  at  sections  i  and  ifl, 
while  forward  differences  are  used  in  place  of  the  derivatives 
of  N,  Q  and  M,.  The  difference  equations  are. 


C  CWT).  -t-CQ  -AQ  =;-DH  -tftNi.  -v2.P, 


C2.^s 

C  C'Ut"),  -vdQ.  -6Q.  =cn  -An.  -V- c  P  . 

c  Jjk  L^k  tk 

C2..f7) 


Sk  ^  i,k  i+\^k  5  ii, 


wV>c\re  C  =  /°b  la. 

‘  T* 

wt"=  £w 


A  ag 


r  S.n  C^_|) 


UT  = 


at^ 
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2.4.4  Futnerical  Method  of  Solution: 

The  solution  is  started  out  from  the  specified 
initial  conditions  and  then  information  is  obtained  at 
successive  time  intervals.  At  apart.  The  method  is  as 
follows: 

(a)  The  arch  is  divided  in  (m-1)  segments  as 

shown  in  Pig.  2v6.  The  ordinates  of  the  radial  and  tangential 

components  of  the  load  pulse  acting  on  each  segment  are  computed 

at  time  intervals  At  apart  so  that  the  values  of  loads  at 

discrete  times  t^,  t^ ,  t^  are  known  as  shown  in 

Fig.  2.10a  or  2.11b.  These  load  values  are  nondimens ion- 

f  13 ) 

alized.  The  natural  period,  T  is  calculated  as  follows  , 


Compression  mode 


1800 


where  T  is  in  seconds 
and  r  is  in  feet. 

O 


Flexural  mode 


•415  000 


-  I  - 


where  T  is  in  seconds 

and  Tj,  and  d  are  in  inches. 
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O  =■  Decaj  ime_  of  ^l^e-  loa4  . 


Figure  2.10  Load-Time  Curves 

(b)  Starting  with  the  undeformed  arch 
equation  (2.62b)  at  time  t^,  the  solution  for  the 
response  is  launched  by  using  the  ’Acceleration-pulse' 
method^ in  this  method  the  following  formulae  are 
Used  for  extrapolating  the  displacements  ^  and 
Ul  k  of  a  secion  1(1  =  2,3> . . .  ,in-l )  at  time  tk 

(k  =  2,3,...). 


2W  -  W.  .-tCWT)  Ca-«-> 


(2-6'^') 
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and  U.  =2U.  -U.  +  cut")  CAt)^ 

i^k-\  L^k-I 


where  At  =  the  nondimens ional  time  interval. 

At  time  tj_ ,  for  any  section  i  (i  =  2,3,..., 
m-1)  the  following  two  types  of  extrapolation  formulae 
are  used. 

(1)  If  the  load-time  curve  is  as 
shown  in  Pig,  2,10a, 


QIn 


V.  =  I  CWt") 

1-1  i'o 


i  u.  =  i  cut':).  caV 

t .1  i  CO 


y  C  2.70  a:) 


J 


(2)  If  the  load-time  curve  is  as 


shown  in  Fig.  S.Jflb, 


-r 


V7.  =  CWT  ).  CAt) 


1,1  ^ 


y  (2. Yob) 


and  U.  =-L  cut").  CAt)^ 

L,1  €  LI 

'  / 
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In  cases  when  equations  (2.70b)  are  used,  sui 
iteration  on  the  values  of  both  (V/T")j^  ^  and 

becomes  necessary. 

(c)  At  time  t^,  equation  (2., 62b) 

indicates  that 

U  ,  =  W.  =  M  .  =  N.  =  Q.  =  O 

i  o  to  t,o 

Using  load  ^  equilibrium  equations  (2.66), 

(2.67)  and  (2,68),  the  values  of  (WT").  and  (UT").  ^ 

J  j  ■) . 

(j  =  li.*.i  -  segments)  are  obtained.  Also, 

approximating  the  radial  and  tangential  accelerations  of 
section  i  by  the  average  of  the  corresponding  acceler¬ 
ations  of  segments  j-1  and  J, 


tf 

It 

(WT  ).  =  -L 

(WT  ). 

+  cwT 

t  k  2_  i 

f  1 

CUT).  -i. 

// 

j 

•t 

CUT  ). 

+  CUT  ). 

k  2- 

J-'.k 

where  (i=  2,3,...,  m-1), 

the  values  of  (WT").  .  and  (UT").  are  obtained. 

1,0  1,0 


C2-VI) 
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(d)  At  time  t^  =  tjj+At, 

( i)  Using  either  equations  (2.70a)  or 
(2,70b)'''  as  need  be,  values  of  ^  and  ^  , 

(i  -  2,...,  m-1)  are  obtained.  The  boundary  conditions 

[^equation  (2.62a)3  require  that. 


W,  ,  =  U,  -  =  M-,  T  =  T  =  U  =  Fl  T  =  0 

1^1  in  f  1  m  ^  JL  rn  ^ 

(ii)  Using  equations  (2.63),  (2,6[j.)  and 

2.65),  values  of  e^^  ^  and  (i  =  l,...,m)  and 

(Adl0).-,  (j  =  l,.,.,  m-1)  are  obtained. 

J  > 

(iii)  For  each  section  i,  equations  (2.55a) 
and  (2.55b)  are  employed  to  obtain  concrete  strains  e^ , 

and  steel  strains  e^ ,  eg. 


'"'when  equations  (2.70b)  are  used,  the  values  of  (WT")^^  ^ 

and  (UT"),  n  (i  =  l,...,m)  are  not  available.  Hence, 

1,1 

certain  reasonable  values  are  assumed;  at  first  operations 
(i)  through  (vii)  are  performed  and  values  obtained  in 
step(vii)  are  compared  with  the  assumed  values.  If  the 
agreement  between  them  is  not  satisfactory,  the  newly 
obtained  values  of  (WT")^  ^  and  (UT")^  are  employed  in 
step  (i)  and  steps  (i)  through.  ( vii)  are  repeated.  This 
process  is  continued  until  a  satisfactory  agreement  between 
the  values  used  in  step  (i)  and  those  obtained  in  step 
(vii)  is  achieved. 
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(iv)  Based  on  the  values  of  ej_  and  e^, 
each  section  is  classified  as  Case  1  [^equations  (2,56), 

(2.57) 3^  Case  2  [^equations  (2.58),  (2.59)3  or  Case  3 
[equations  (2.60),  (2. 6^)3.  Also,  strains  e^  and  e^ 

for  each  section  i,  are  conpared.  with  the  yield  strain  of 

steel  e^y  and  are  replaced  by  e^^  (with  proper  sign) 

if  found  to  be  greater  than  e^^  in  magnitude. 

(v)  Making  use  of  the  equations  (2.56), 

(2.57)  or  (2.58) ,  (2.59)  or  (2.60),  (2.61)  values  of 
S'i  ^  and  ^  ( i  =l,...,ra)  are  obtained. 

(vi)  Loads  P  .  ,,  P  .  ,  and  (  AdO) 

rj  e  J  >-1-  J 

calculated  in  step  (ii),  (j  =1 , . . . ,m-l )  and  equations 
(2.66),  (2.67),  (2.68)  are  used  to  calculate  (WT"),-  n, 

J  >  i- 

(UT")j  1  (j  =  l,...,m-l)  and  ^  (i  =  l,...,m).  For 
the  (m-l)  segments,  (3™-‘3)  equations  are  now  available 
while  the  unknowns  form  a  total  of  (3m-2),  i.e.,  (m)  Qs, 
(m-l)  WT"s  and  (m-l)  UT"s.  In  order  to  overcome  this 
difficulty,  it  is  assumed  that 

(WT")j  1  9  where  j  =  1  and  m-l 

J  ^ 

This  assumption  implies  that  the  radial  accelerations  of 
segments  nearest  to  the  two  supports  are.  neglected.  This 
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leaves  quantities  as  unknown.  Excluding  equation  (2,66) 

as  applied  to  segment  (m-l),  a  total  of  {3ra-4)  equations  are 

available.  Solving  these  equations,  values  of  (WT").  ,, 

J  ^  -L 

(UT")j  "  2,,..,m-2)  and  ^  ^  ~  are  obtained. 

However,  in  the  case  of  uniformly  distributed  symme¬ 
trical  loading,  making  use  of  the  fact  that  the  shear  Q  at 
the  crown  is  zero,  equations  (2,66),  (2.67)  and  (2,68)  are 
simultaneously  solved  for  half  the  arch.  The  solutions  of 
of  these  equations  yield  the  values  of  (WT"),  (UT")  and  Q 

for  one-half  of  the  arch:  the  corresponding  values  for  the 
other  half  are  obtained  by  using  symmetry, 

(vii)  Finally,  using  equations  (2.71)  for  k  =  1, 

values  of  (WT")j_  and  (UT")j^  1  ~  2 , .  .  .  ,m-l-seotions )  are 

obtained. 

At  the  end  of  step  (vii)  all  the  information  regarding 
the  internal  forces  M,  N  and  Q,  the  radial  and  tangential  dis¬ 
placements  W  and  U,  and  the  corresponding  accelerations  WT" 

and  UT"  is  available  at  time  t  for  all  the  m  sections. 

1 

(e)  At  each  successive  time  t  ,  t  tv  the  seven 

steps  of  (d)  are  used  to  obtain  all  the  information  about  the 
arch  at  the  discrete  sections,  the  only  difference  being  that 
in  step  ( i)  equations  (2,69)  are  used  to  calculate  displace¬ 
ments  Wj^  ^  and  U^  (i  =  2,,,.,  m-l  -  sections). 

The  steps  (a)  through  (e)  thus  give  the  response 
of  the  arch  to  a  time-dependent  load. 

2,4.5  Calculation  of  Ultimate  Load: 

The  ultimate  load  of  an  arch  is  once  again 
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leaves  (3tii-4)  quantities  as  unknown.  Excluding  equation  (2.66) 

as  applied  to  segment  ( m-1 ) ,  a  total  of  equations  are 

available.  Solving  these  equations,  values  of  (WT")  .  , 

J  >  -*- 

(UT").  (j  "  2,...,m-2)  and  Q,  ,  (1  =  l,...,m)  are  obtained, 

J  >  -L  1,1 

However,  in  the  case  of  uniformly  distributed  symme¬ 
trical  loading,  making  use  of  the  fact  that  the  shear  Q  at 
the  crown  is  zero,  equations  (2.66),  (2,67)  and  (2.68)  are 
simultaneously  solved  for  half  the  arch.  The  solutions  of 
of  these  equations  yield  the  values  of  (WT"),  (HT")  and  Q 

for  one-half  of  the  arch:  tti  e  corresponding  values  for  the 
other  half  are  obtained  by  using  symmetry, 

(vli)  Finally,  using  equations  (2,71)  for  k  =  1, 

values  of  (WT")j_  ^  and  (UT")j^  ^  (i  =  2, . .  ,  ,m-l-sections )  are 

obtained . 

At  the  end  of  step  (vii)  all  the  information  regarding 
the  internal  forces  M,  N  and  Q,  the  radial  and  tangential  dis¬ 
placements  W  and  U,  and  the  corresponding  accelerations  WT" 
and  UT"  is  available  at  time  t  for  all  the  m  sections. 

X 

(e)  At  each  successive  time  t  ,  t  tv  the  seven 

3  3 

steps  of  (d)  are  used  to  obtain  all  the  information  about  the 
arch  at  the  discrete  sections,  the  only  difference  being  that 
in  step  (i)  equations  (2.69)  are  used  to  calculate  displace¬ 
ments  ^  and  U^  (i  =  2,,..,  m-1  -  sections). 

The  steps  (a)  through  (e)  thus  give  the  response 
of  the  arch  to  a  time-dependent  load. 

2,4.8  Calculation  of  Ultimate  Load: 

The  ultimate  load  of  an  arch  is  once  again 


defined  to  be  that  load  under  which  the  maximum  compressive 
strain  in  concrete  exceeds  the  ultimate  strain  e^.  There¬ 
fore,  as  explained  in  §  2.3.5,  this  method  gives  ultimate 
loads  for  certain  distributions  of  loading  such  as 
(a)  uniformly  distributed  symmetrical  loads,  (b)  uniformly 
distributed  antisymmetrical  loads,  and  (c)  symmetrical  and 
antisymmetrical  concentrated  loads,  except  a  concentrated 
load  at  the  crown. 

For  a  given  distribution  of  loading  and  for  a 
given  load  pulse  such  as  the  one  shown  in  Pig.  2.10b 
(defining  t  ,  t„  and  t,)  the  ultimate  load  is  obtained 
as  follows; 

(1)  Starting  at  a  low  value  of  the  peak  load 
p  the  arch  is  analysed  at  discrete  times  t^,  tj_,...,tj^ 
using  the  procedure  outlined  in  §2,4.4.  At  each  time  the 

maximum  v'alue  of  the  concrete  compressive  strain  is  compared 
with  the  concrete  ultimate  strain  to  check  for  failure. 

The  analysis  is  continued  either  until  t^  =  t^  (if  t^  T) 
or  until  -  T  ( if”  <  T) .  Since  the  maximum  response 

of  the  arch  occurs  at  a  time  t^  <  T  (usually  0.5t  to 

0.75T)^^^^,  the  analysis  of  the  arch  upto  T  is 

sufficient  for  investigation  of  the  possibility  of  failure 
under  the  load  p  . 
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(2)  The  peak  load  p  la  increased  by 

t 

a  certain  percentage  and  the  analysis  as  explained  in 
step  (1)  is  repeated. 

(3)  This  process  of  increasing  the 

peak  load  p^  and  analysing  the  arch  is  continued  until 
the  ultimate  strain  e  is  exceeded.  The  value  of  p 
at  which  this  excessive  strain  is  produced  is  the  ultimate 
load  of  the  arch  for  a  given  distribution  of  loading  and 
a  given  load-time  dependence. 

2.4.. 6  Digital  Computer  Program: 

A  digital  computer  program  prepared  to  perform 
the  calculations  involved  in  the  methods  outlined  in 
articles  2.I4..I4-  and  2.[|..5  ^3  presented  in  Appendix  IV. 

A  flow-chart  of,  the  program  is  also  presented  in 
Appendix  IV. 

2.1|.*7  Selection  of  Space  and  Time  Intervals: 

Before  the  numerical  method  presented  in 
2. 4.. 4  can  be  used,  it  is  necessary  to  determine  the 
approximate  values  for  the  parameters  a©  and  A.t  used 
for  discretization  of  the  continuous  system.  As  ln§2.3.7» 
again  a  semicircular  arch  is  divided  into  20  segments  so 
that  A0  is  equal  to  nine  degrees. 
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In  order  to  establish  the  time  interval  At, 
information  regarding  the  problem  of  wave  propogation  in 
a  one  dimensional  elastic  system  proves  to  be  useful. 
Crandall^ has  shown  that  in  a  one  dimensional  elastic 
system,  the  time  interval  At  is  related  to  the  ratio  of 
the  space  interval  AX  and  the  seismic  velocity,  in 
the  medium.  This  relation  says  that  for  the  stability 
of  the  numerical  method  of  solution. 


At 


Ci-ya.) 


Crandall  has  also  suggested  that  if  At  is  much  smaller 
A  X 

than  p —  ,  the  accuracy  of  the  numerical  method 

deteriorates . 


Wivv, 


wy>cl 


AX  =  r  A0 

O 


> 


(2.73) 


Semicircular  arches  of  Type  A  (refer  to  §  3*3)  under 
different  distributions  of  loading  with  different  load¬ 
time  functions  have  been  analysed.  It  is  found  (Pig.  2.11) 
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CM 

O 

d 


o 

d 


that  for  At 


the  convergence  is 


=  AX  or  ^  ^ 

adequate  and  satisfactory  deflection-time  curves  are 

Ax 

obtained.  However,  for  At  -g^the  deflection-time 
curves  show  erroneous  trends.  Therefore,  for  all  the 

A  V 

dynamic  cases  the  criterion  At  = is  used#  permitting 

'"s 

the  use  of  the  largest  time  interval  without  causing  numer¬ 
ical  instability. 

2.5  CONVENTIONAL  THEORY  FOR  PREDICTING  ULTIMATE  LOADS 
AND  DEFLECTIONS- 


In  this  section  a  brief  review  is  given  of  the  con¬ 
ventional  methods  of  predicting  the  ultimate  loads  of 
arches  under  (1)  uniformly  distributed  symmetric  radial 
load,  (2)  concentrated  load  at  the  crown  and  (3)  anti¬ 
symmetric  loads  at  the  quarter  points.  The  expressions 
for  deflection  at  the  points  where  the  deflections  are 
measured  in  tests  (refer  section  3*l4-*2  (c)  in  following 
chapter),  are  also  given.  The  Initial  tangent  modulus  of 
concrete  is  to  be  used  in  these  expressions.  The  expressions 
for  deflections  are  derived  by  assuming  linearly  elastic 
material,  using  the  usual  virtual  work  approach. 

2.5.1  Uniformly  Distributed  Symmetric  Radial  Load: 

(a)  Ultimate  Load:  If  the  secondary  bending 
effects  are  neglected,  the  ultimate  load  is  given  by  the 
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following  expression. 


u 


TH 


(0-85  I  -V  f  M 

c  ^  ^ 


O 


(2.yA) 


where 


t  =  Total  depth  of  arch  in  inches 

•p 

°  =  Radius  of  arch  in  inches 

b  =  Width  of  arch  in  inches 

p  =  Ultimate  load  in  Ibs/lnch 

^TH  A„  +  A  ' 

3  3 

p,  =  Total  steel  percentage  =  - 

^  bt 


(b)  Deflection:  The  radial  deflections  of  the  arch 
at  crown  and  at  points  5i+°  from  the  supports  are  given  by 
the  following  expressions. 


W 

C 


w 


\.S^6  TH  0 
AE. 


\-10 


r 

TH  O 

AE. 


CZ.T5) 


C2-76) 


where 


A  =  A^  [1  +  (n-1)  p^] 

E  =  (1800000  +  390  f)  in  bSL. 

c  c 

^TH  ~  Distributed  load  in  Ibs/inch 

A  =  Area  of  concrete  section  in  sq.  inches 

c 

E 

n  =  Modular  ratio  =  ^  =  10  (assumed) 

^c 

w  --Radial  deflection  at  the  crown  in  inches 

c 

w  =  Radial  deflection  at  $h.°  points 
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2.5-2  Concentrated  Load  at  Crown: 

(a)  Ultimate  Eoad:  The  ultimate  load  of  the 
arch  under  concentrated  load  can  be  calculated  by  using  the 
conventional  plastic  theory  as  the  test  specimens  of  the 
test  program  are  underreinforced.  It  can  be  shown  that 
the  arch  will  collapse  with  the  formation  of  three  hinges, 
one  at  the  crown  and  two  others  at  38°  from  the  supports. 

It  is  not  easy  to  derive  an  explicit  formula  for  the 
ultimate  load  because  the  effsob  of  thrust  on  the  moment 
capacity  of  the  section  has  to  be  considered.  An  iterative 
procedure  described  below  can  be  used. 

In  the  first  cycle,  the  ultimate  load  can 
be  calculated  by  using  equilibrium  conditions  and  neglecting 
the  effect  of  thrust.  Two  unknowns  -  horizontal  reaction  H 
and  ultimate  load  can  be  found  by  using  the  conditions 
that  the  moments  at  two  hinges  should  be  the  ultimate  moment 
M^,  which  is  given  by 

M  =  A  ?  Cd  -  fl/zD  Cz.Tv) 

u  S  y 


wWtK 


(o-85  f'  b) 


Once  this  ultimate  load  is  calculated,  it  can  be  used  to 
calculate  the  thrust  at  the  crown  and  at  the  38°  points. 


With  these  values  of  thrusts,  new  values  of  moment  capacities 
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at  crown  and  at  38°  points  can  be  found  from  the  graph 


of  vs.  M/M^.  (See  for  example  Pig.  Sk  -  10.2, 

Reference  No.  13<  in  bibliography.)  With  these  values  of 
moment  capacities  a  second  cycle  of  iteration  will  give 
a  new  value  of  H  and  P^.  This  process  can  be  repeated  until 
a  satisfactory  convergence  is  obtained.  Usually  two  or 
three  cycles  should  give  satisfactory  results. 

(b)  Deflection:  The  expression  for  the 
radial  deflection  at  crown  is  as  follows. 

P 

w  =  O.OO-I  ri‘^8) 

where 

E 

c 

p 

TH 
w 

c 

moment  of  inertia  of  reinforced  concrete  section  should  be 
used.  However,  as  recommended  in  references  11  (Section  2.6) 
and  19  in  bibliography,  it  is  decided  to  use  an  average 
of  the  moments  of  inertia  of  the  cracked  and  uncracked 
sections  (l.e.,  • 


c 

=  Average  of  the  moments  of  inertia  of 

the  cracked  and  uncracked  sections  in  in.^ 

=  (1800000  +  390  f)  in  psi 

c 

=  Concentrated  load  at  crown  in  pounds. 

=  Radial  deflection  at  crown  in  inches. 

It  is  uncertain  as  to  what  value  of  the 


I 
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2.5.3  Antisymmetric  Concentrated  Loads  at  Quarter 
Points : 

(a)  Ultimate  Load:  The  ultimate  load  in  this 
case  can  be  calculated  also  by  using  plastic  theory.  Under 
this  type  of  loading  it  can  be  shown  that  the  arch  collapses 
by  the  formation  of  two  hinges,  one  at  each  quarter  point. 

To  take  into  account  the  effect  of  thrust  on  the  moment 
capacity  of  the  section,  the  following  iterative  procedure 
can  be  used. 

In  the  first  cycle  of  iteration,  if 
the  effect  of  thrust  is  neglected,  the  ultimate  load 
is  given  by 


P 

u 


2  M 

U 

)r 


C2.Y‘^) 


where  =  A  ^  Cd  -  Q/z) 

u  ^  J 

and  (X  —  - - — 

0.85  f'b 

C 

With  this  value  of  ultimate  load,  thrust 
at  two  quarter  points  can  be  calculated.  These  thrust 
values  can  be  used  to  determine  the  new  lalues  of  moment 
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deflection  is  calculated  by  using  formulae  presented 
in  §  2.5  and  the  analytically  predicted  ultimate  load. 
The  experimental  value  of  is  calculated  by  using  the 
value  of  the  maximum  deflection  under  the  ultimate  load 
measured  as  described  in  §  3.4»2  (c)  and  the  elastic 
deflection  obtained  by  using  again  the  formulae  given  in 
J  2.5  and  the  experimentally  determined  ultimate  load. 

2.7  COMPARISON  OP  DYNAMIC  ULTIMATE  LOADS  -  NONLINEAR 
THEORY  AND  APPROXIMATE  THEORETICAL  ANALYSIS 
The  nonlinear  theory  developed  so  far  is  used  to 
obtain  ultimate  loads  of  arches  under  a  triangular 
dynamic  load  pulse  with  zero  I’ise  time  (Pig.  2.12). 

The  distributions  of  loading  along  the  arch  are 

a)  Compression  mode  type 

and  b)  Deflection  mode  type 

The  geometry  and  material  properties  of  the  arches 
is  as  follows: 

b  =  10"  ,  t  =  15!"  r^  =  300" 

Jlf  ^  =  90°  ,  p^  =  0.025  d  =  13.5" 

f •  =  4000  psi,  f ,  =  40,000  psi, 

dc  ^ 

E„  =  30  X  10^  psi 

3 
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Figure  2.12  Triangular  Load  Pulse 

Various  values  of  the  ultimate  loads,  are  obtained 
for  different  duration  of  loading  t|,,  with  both  types 

of  distributions. 

Approximate  analysis  based  on  a  single-degree 
freedom  system  is  also  conducted  for  the  same  arches. 
The  following  formulae^^^  ^  ^  (19)  used  for  this 


purpose: 


Compression  mode  loading 


Pm  ~ 


-L.,  2.^.1, 

Ati 


where  Pj„  =  Dynamic  load  (ibs/ln) 


r  =  (0.85  f’  Pt  fdy)  M  (Ibs/ln) 


Natural,  period  of  the  arch  in  compression 
mode  (milliseconds) 


21.6 

Mean  radius  of  the  arch  (Inches) 
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Ductility  factor 


F  = 


Here,  is  ca'lculated  for  two  values  of  fx,  viz. 
p.  =1.3  and  2,  and  various  values  of  . 


T 

-  (2 

A  ti 


T  =  Natxiral  period  of  the  arch  in  flexural 
mode  (milliseconds) 

(V  +  1.5 
(  V  00  )^-  1  - 

Here  various  values  of  t^  and  two  values  of 
vix.  /X  =  2  and  3  are  used  for  calculating  p^. 

The  analytical  results  obtained  from  the  nonlinear 
theory  are  compared  with  results  obtained  from  the  approx¬ 
imate  analysis  (20).  This  comparison  is  shown  in  Pig.  2.13. 

There  is  a  good  agreement  between  these  results  when  the 
duration  of  loading  t^  is  greater  than  about  half  the 


(  ^0  •  ^o)^ 
425  d  sj~v 
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natural  period,  T.  However,  the  approximate  analysis 
seems  to  give  unconservative  results  for  t^^  less  than 
about  0.5  T.  A  possible  explanation  for  such  results 
may  lie  in  the  fact  that  under  loads  of  short  duration, 
more  than  one  mode  of  vibration  are  excited  and  analysis 
based  on  a  single-degree  freedom  system  becomes  inadequate. 
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FAILURE  PEAK  LOAD  -  (LBS/INCH)  NONDIMENSIONAL  FAILURE  PEAK  LX)AD- 


i 


0  0.5  1.0  1.5  2.0 

NONDIMENSIONAL  DURATION  OF  LOAD  - 


(a)  COMPRESSION  MODE  LOADING 


0  0.5  1.0  1.5  2.0  . 

NONDIMENSIONAL  DURATION  OF  LOAD--^ 


(b)  DEFLECTION  MODE  LOADING 

FIGURE  2.13- COMPARISON  OF  ANALYTICAL  RESULTS  WITH 
APPROXIMATE  ANALYSIS 
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CHAPTER  3 


EXPERIMENTAL  INVESTIGATION 
3.1  TEST  SPECIMENS 

All  test  speciraens  are  reinforced  concrete,  semi¬ 
circular  arches  of  radius  18"  and  of  rectangular  cross- 
section.  For  the  purpose  of  these  tests,-  two  sizes,  one 
with  a  cross-section  of  2"  x  2"  and  the  other  with  a 
cross-section  of  1”  x  2"  were  chosen.  The  purpose  of 
choosing  two  sizes  was  to  detect  if  possible, the  in¬ 
fluence  of  the  natural  period  on  the  dynamic  response  of 
the  arch  for  a  given  rise  time  of  the  dynamic  load.  The 
details  of  the  geometry,  cross-section  and  reinforoement 
of  the  arches  are  shown  in  Table  3*1  and  Figs.  3.I  and  3.2. 
The  specimens  with  2"  x  2"  cross-section  are  designated 
as  Type  A,  and  those  with  1”  x  2”,  as  Type  B. 


TABLE  3*1 


Specimen 

Type 

Cross- 

section 

Reinforcement 

A 

3 

(  sq . in. ) 

A 

2"  deep  x 
2"  wide 

l|  #7  wires 
(diameter  =  0.177”) 

0. 0492 

0.01+92 

B 

1"  deep  X 

2  wide 

[|.  #12  wires 
(diameter  =  O.IO^” 

0.0173 

0.0173 
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Suitable  wooden  forma  were  used  to  cast  the  speciraena. (Pig. 3. 3) 
A  aimple  device  aa  ahown  in  Fig.  3*4  waa  uaed  to  make  the 
reinforcement  cagea.  The  wirea  were  firat  bent  to  the 
required  radiua  in  a  wire  bending  machine^  and  then 
inaerted  into  the  device  mentioned  above.  Stlrrupa  were 
tied  at  predetermined  apacing. 

The  concrete  waa  mixed  for  2  to  3  minutea  in  a  nine 
cubic  feet  capacity  tilting  drum  type  mixer.  The  apecimena 
were  left  in  the  formwork  for  about  48  houra,  after  which 
they  were  removed  from  the  forma  and  cured  in  the  air  of 
the  laboratory  until  teated.  With  each  apecimen,  three 
6"  X  12"  control  cyllndera  were  alao  caat  and  cured  under 
the  aame  conditlona.  Theae  control  cyllndera  were  uaed 

I 

to  determine  the  ultimate  compreasive  atrength  f^  of  the 
concrete,  at  the  time  the  apecimen  waa  teated. 

3.2  MATERIAL  PROPERTIES 
3*2.1  Concrete: 

I 

The  concrete  mix  waa  dealgned  for  an  f  ^  of 

3000  pai  at  aeven  daya.  The  proportiona  by  weight  of  the 

cement,  aand,  gravel  and  water  were  aa  followa: 

1  part  by  weight  of  high  early  portland 
cement# 

2.64  part  by  weight  of  aand, 

2.10  part  by  weight  of  coarae  aggregate, 

7.25  giillona  of  water  per  aack  of  cement. 
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with  these  proportions  it  was  expected  to  obtain 

equal  to  3000  psi  at  seven  days.  However,  for  some  test 

specimens,  the  actual  f’  obtained  from  the  test  cylinders 

c 

was  higher  because  it  was  not  possible  to  test  the  specimens 
exactly  at  seven  days.  In  some  cases,  specimens  were  tested 
after  as  much  as  thirty  days,  due  to  unexpected  delays.  The 
actual  value  of  f^  of  each  specimen  is  given  in  Chapter  4- 
The  sand  used  had  a  fineness  modulus  of  about  2.5. 
The  maximum  aggregate  size  for  Type  A  specimens  was  3/8", 
while  for  Type  B  specimens  3/l6". 

3.2.2  Steel: 

Black  annealed  mild  steel  wires  were  used  as 
the  reinforcement.  For  Type  A  specimens,  #7  wires  (diameter 

0,177")  were  used,  while  for  Type  B  Specimens  #12  wires 
(diameter  =  0.105")  were  used.  The  strength  properties  of 
these  wires  were  as  follows: 

(a)  #7  Wires: 

Static  ylelu  stress  f^  =  4P  ksi 
Static  ultimate  stress  f-u  =  5l  ksi 
Modulus  of  elasticity  E_  ■=  34-1  x  10®  ksi 

O 

The  stress-strain  curve  of  this  wire  is  given 


in  Fig.  3.5. 


(b)  #12  Wires; 

Static  yield  stress  fy  =  40  l^si 

Static  ultimate  stress  f-u  =  49.5  ksi 

Modulus  of  elasticity  E„  =  31.6  x  10®  ksi 

s 

The  stress-strain  curve  of  this  wire  is  given 
in  Fig.  3.6.  • 

3.3  LOADING  CONDITIONS 

In  this  test  program,  static  and  dynamic  tests  were 
conducted  for  the  following  three  types  of  loading  conditions. 

Type  I  -  Uniformly  distributed  symmetric  radial  load; 

This  was  simulated  by  ten  point  loads, 
spaced  at  equal  intervals  on  the  periphery  of 
the  arch  as  shown  in  Fig,  ,3.7.  Each  interval 
subtends  an  angle  of  l8°  at  the  center. 

Type  II  -  Concentrated  load  at  the  crown  of  the  arch. 
Type  III-  Antisymmetric  concentrated  loads  at  quarter 
points  on  the  arch. 

These  three  loading  conditions  are  shown  in  Fig.  3.7. 

3.4  EXPERIMENTAL  SET  UP 

3,4.1  Modification  of  the  Existing  Loading  Machine; 

/  D  \ 

The  dynamic  loading  machine  ,  constructed 
by  the  Department  of  Civil  Engineering,  M.I.T.,  under  the 
contract  DA-49-129 -ENG-325,  had  to  be  modified  to  suit 
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the  loading  conditions  described  in  section  3*3«  With  the 

earlier  version  of  the  machine  it  was  possible  only  to  apply 

a  single  point  load  on  the  test  specimen.  The  modification 

essentially  consisted  of  mounting  ten  jacks  on  a  suitable 

frame,  and  in  connecting  these  jacks  through  rubber  hoses 

to  the  auxiliary  oil  reservoirs.  In  what  follows,  a  short 

description  of  loading  jacks,  rubber  hoses,  oil  reservoirs 

and  the  details  of  their  connection  with  the  previously 

existing  apparatus  is  given.  For  a  detailed  description  of 

the  components  of  the  original  machine  and  the  principal 

fl8) 

of  its  operation (reference  la  made  to 

(a)  Loading  Jacks:  The  preliminary  estimate  of 
the  dynamic  resistance  of  the  Type  A  arch  specimens  showed 
that  the  failure  load  under  compression  mode  loading  (i.e.. 
Type  I  loading  condition)  would  be  about  6.5  to  7  kips  per 
jack.  Therefore,  it  was  decided  to  design  the  loading  jacks 
to  develop  a  load  of  12  kips  at  the  maximum  working  oil 
pressure  of  I4.OOO  psi.  The  details  of  design  were  based 
on  considerations  similar  to  those  given  in  reference  (18) 
(Chap.  V,  Art  c).  It  was  found  that  Hannifin  Series  "H” 
square  type  hydraulic  cylinders  (manufactured  by  Haniffln 
Company,  Des  Plaines,  Illinois)  with  the  following  re¬ 
quirements,  were  suitable  for  the  purpose  and  accordingly 
adopted . 
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2.  Piston  rod  diameter 


II 


3-  Stroke  =  Ij." 

I4..  Maximum  Working  Pressure  =  I4.OOO  psi 
These  jacks  are  double  acting  -  double  rod  end,  with  cushions 
on  both  ends  of  the  jack. 

(b)  Reservoir:  In  the  modified  loading  system, 
two  reservoirs  are  used.  One  of  them  is  interposed  between 
the  pump  and  the  push  side  of  the  loading  jacks  while  the 
other  is  placed  between  the  pull  side  of  the  jacks  and  the 
dump  valve  of  the  original  dynamic  loading  machine.  These 
reservoirs  are  in  the  form  of  cylindrical  containers  with 

a  sufficient  number  of  outlets  which  can  be  connected  to 
the  jacks.  These  reservoirs  are  shown  in  Pigs.  3.8  and 
3.9.  Various  details,  such  as  the  diameter  of  cylindrical 
container,  number  of  outlets,  etc.,  are  also  marked  in 
Pigs.  3*8  and  3*9*  The  reservoirs  are  designed  to  with¬ 
stand  an  oil  pressure  of  5000  psi.  The  main  purpose  of 
these  reservoirs  is  to  provide  the  means  of  supplying  oil 
to  all  the  jacks  at  the  same  pressure  so  that  the  Jacks 
can  develop  equal  loads. 

(c)  Rubber  Hoses:  In  order  to  have  the  flexibility 
in  connections  between  the  reservoirs  and  the  jacks,  rubber 
hoses  were  used  instead  of  metal  tubing.  These  hoses  are 

3  wire  braid  high  pressure  hoses  (working  pressure  =  5000  psi. 


1" 

and  bursting  pressure  =  20,000  psi)  with  ^  internal  diameter. 

(d)  Hydraulic  Connection:  The  hydraulic  connec¬ 
tions  between  various  components  of  the  loading  system 
are  shown  in  Fig.  3*10. 

The  main  cylinder  of  the  machine  is  connec¬ 
ted  to  the  push  side  reservoir  which  in  turn  is  connected 
to  the  push  side  of  the  jacks..  The  main  dump  valve  is 
connected  to  the  pull  side  reservoir  which  in  turn  is 
connected  to  the  pull  side  of  the  jacks.  Heavy  duty  steel 
pipes  are  used  for  connecting  the  main  cylinder  with  the 
push  side  reservoir  and  the  dump  valve  with  the  pull  side 
reservoir.  Connections  between  the  reservoirs  and  the 
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sections.  Fig.  3.11  diows  this  frame  with  ten  jacks 
mounted  for  Type  I  loading.  For  Type  II  loading ^ only 
one  jack  is  mounted  at  the  proper  place  on  the  frame  as 
shown  in  Fig.  3.12.  Two  jacks  are  mounted  on  the  frame  as 
shown  in  Fig,  3.13  for  Type  III  loading.  In  order  to  ensure 
that  the  applied  load  remains  radial,  as  the  arch  deforms 
under  load,  the  jacks  are  mounted  so  as  to  allow  them  to 
rotate  with  the  arch.  This  is  done  by  using  pin  connections 
between  the  jacks  and  the  frame  and  between  jacks  and  the 
specimen. 

The  pinned  supports  for  the  test  specimen  is 
simulated  by  using  specially  fabricated  devices  as  shown 
in  Figs.  3.14,  3.15,  3.16  and  3.17. 

3.4*2  Instrumentation  and  Measurements: 

In  order  to  study  the  behavior  of  the  test 
specimens,  it  is  necessary  to  measure  the  following 
quantities : 

(a)  Applied  Loads, 

(b)  Reactions, 

(c)  Radial  deflection  at  crown  and  at  2,  points, 

4 

(d)  Natural  frequency  of  the  test  specimens. 
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(a)  Applied  Loads;  The  applied  loads  were 
tneasiired  by  mounting  suitable  load  cells  on  the  piston  rod 
of  the  jacks.  The  range  of  failure  loads  between  type  A 
specimens  under  compression  mode  loading  (Type  I)  and 
type  B  specimens  under  deflection  mode  loading  (Type  III) 
was  quite  large.  Type  A  specimen  under  compression  mode 
had  an  estimated  failure  load  of  about  7000  lbs.  per  jack, 
while  the  estimated  failure  load  for  type  B  specimen  under 
deflection  mode  was  as  low  as  100  lbs.  per  jack.  Therefore 
three  different  types  of  load  cells  had  to  be  made  for  each 
of  the  three  load  types.  These  load  cells  are  shown  in 
Pig.  3.18. 

Load  cells  for  Type  I  and  Type  II  loadings  were 

essentially  hollow  circular  aluminum  rods  on  which  are 

mounted  eight  C-7  strain  gages  (4  active  gages  and  4  dummy 

gages,  resistance  of  each  gage  =  500  ohms).  Load  cells 
for  Type  III  loading  were  in  the  form  of  a  u-shaped 

aluminum  frame  on  the  vertical  sides  of  which  are  mounted 

eight  C-7  strain  gages.  These  strain  gages  are  connected 

to  form  a  suitable  Wheatstone  bridge  circuit.  The  signal 

from  this  strain  gage  bridge  was  fed  into  an  eighteen 

channel  recording  oscillograph.  Type  5-ll4“P3,  manufactured 

by  the  Consolidated  Electrodynamics  Company.  (Hereafter 

in  this  report  this  equipment  will  be  referred  to  as  the 

C.  E.  C.  recorder).  The  traces  of  the  galvanometers 
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of  the  C.E.G.  recorder  were  recorded  on  photographic  paper 
in  both  dynamic  and  static  tests. 

Initially  a  few  pilot  tests  were  made  on 
type  ''a"  specimens  under  Type  I  loading  and  in  these  tests 
all  ten  jack  loads  were  measured  to  determine  whether  the 
loads  developed  by  all  Jacks  were  equal  or  not.  These 
pilot  tests  did  confirm  that  the  loads  were  equal  and  as 
a  consequence,  it  was  decided  to  measure  only  five  jack 
loads  during  the  tests. 

(b)  Reactions:  Vertical  and  horizontal  reactions 
were  measured  by  the  load  cells  as  shown  iti  Pig»  Again 

for  the  reasons  mentioned  above  [article  3*l4.*2  (a)]  various 
types  of  load  cells  had  to  be  made.  These  load  cells  are 
in  the  form  of  solid  or  hollow  aluminum  rods  of  suitable 
diameter,  with  square  aluminum  plates  at  their  ends.  These 
load  cells  were  useful  for  measuring  relatively  large  re¬ 
actions.  The  ring  type  load  cells  which  essentially  con¬ 
sisted  of  an  aluminum  ring  between  two  square  plates  were 
used  to  measure  relatively  small  reactions.  Eight  C-7 
strain  gages  were  mounted  at  suitable  places  on  these  load 
cells  as  shown  in  Pig.  3*19.  The  gages  were  connected 
to  form  a  suitable  Wheatstone  bridge  circuit.  The  output 
of  the  load  cell  was  measured  by  the  C.E.p.  recorder  in 
both  static  and  dynamic  tests. 


-81- 


The  power  to  all  the  load  cells  was 
supplied  by  a  Sorenson  Transistorized  D.C.  Power  Supply 
(Model  No.  QR-36-UA).  The  range  of  output  voltage  and 
output  current  of  this  power  supply  unit  is  0  -  36  volts 
and  0-4  amps,  respectively,  with  a  ia  regulation  of 
0.02. 

(c)  Measurement  of  Deflections:  In  all  tests, 
radial  deflections  of  the  test  specimen  were  measured  at 
the  following  points. 

For  Type  I  loading,  radial  deflections 
were  measured  at  the  crown  and  at  the  points  54.  degrees 
from  the  supports  (Fig.  3.20).  Deflections  were  measured 
at  these  54'^  points  rather  than  at  the  quarter  points 
because  it  was  difficult  to  attach  the  deflection  measuring 
device  at  the  latter  due  to  the  presence  of  the  loading 
jack.  In  Type  I  loading,  it  was  found  during  pilot  tests 
that  the  supporting  frame  also  deflected  appreciably  and 
therefore  deflections  of  the  support  were  also  measured. 
Electric  inductance  gages  of  the  moving  core  solenoid 
type,  commonly  known  as  Linear  Variable  Differential 
Transformers  (L.V.D.T.),  were  used  to  measure  the  deflections. 
These  gages  were  of  the  type  1000  -  SS  -  L  (manufactured 
by  Schaevitz  Engineering  Corportatlon)  with  a  linear  range 
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of  +  1.0".  The  core  and  the  transformer  of  the  L.V.D.T. 
were  mounted  on  a  special  attachment  shown  in  Fig.  3 >21. 

One  end  of  this  attachment  was  connected  to  the  specimen, 
while  the  other  end  was  connected  at  a  suitable  point 
independent  of  the  supporting  frame.  These  connections  were 
such  as  to  allow  the  attachment  to  rotate  so  that  it  re¬ 
mained  radial  when  the  arch  deformed  under  load.  Pig*  3*11 
shows  the  set  up  of  L.V.D.T. 's. 

For  Type  II  loading,  radial  deflections 
were  measured  at  the  crown  and  at  two  quarter  points  (Fig.  3*20). 
At  quarter  points  the  attachments  as  described  in  the 
above  paragraph  were  used  to  mount  the  L.V.D.T. 's  while 
at  the  crown  the  L.V.D.T.  was  mounted  directly  on  the 
jack.  The  setup  of  L.V.D.T. 's  is  shown  in  Pig.  3.12. 

For  Type  III  loading,  radial  and  tangential 
deflections  were  measured  at  the  two  quarter  points.  Type 
2000  -  S3  -  L  L.V.D.T. 's  with  a  linear  range  of  +  2.0" 
were  used  to  measure  the  radial  deflect  long,  while  tangential 
deflections  were  measured  by  type  1000  -  SS  “  L  L.V.D.T. 's. 

The  L.V.D.T. '3  measuring  the  radial  deflections  were  mounted 
directly  on  the  jack  while  the  L.V.D.T. 's,  measuring  tan¬ 
gential  deflections,  were  mounted  on  an  attachment  as 
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shown  in  Fig.  3.13. 

(d)  Measiirement  of  Natural  Period:  The 
natural  periods  of  types  A  and  B  specimens  in  the  first 
bending  mode  (i.e„,  antisymmetric  mode)  and  in  the  second 
bending  mode  (i.e.,  mode  corresponding  to  the  configuration 
of  the  arch  under  a  point  load  at  the  crown)  were  measured. 
For  this  purpose  displacements  corresponding  to  each  mode 
were  given  to  the  arch  specimens,  released  suddenly,  and 
the  resulting  vibrations  measured  by  very  sensitive  L.V.D.T,* 
(Type  020  MS  -  L,  linear  range  -  ±  0,02")  connected  to  the 
C.E.C.  recorder. 

3.5  TESTING  TECHNIQUE 

In  this  section  a  brief  description  of  the  method  of 
static  and  dynamic  testing  is  given. 

3.5.1  Static  Tests: 

The  specimens  were  mounted  in  the  frame  and  con¬ 
nected  to  the  loading  jacks  as  shown  in  Figs,  3.11)  3.12 
and  3,13.  The  specimens  were  loaded  by  continuously 
building  up  the  oil  pressure  on  one  side  of  the  jacks,  until 
the  specimens  failed.  A  continuous  record  of  all  the 
measurements  was  obtained  by  running  the  C.E.C.  recorder 
at  the  slow  speed  of  0,80  in/sec.  The  duration  of  test 
on  an  average  was  3  to  4  minutes. 


3.5*2  Dynamic  Testa: 

The  specimens  were  mounted  in  the  frame  and 
connected  to  the  loading  jacks  as  in  static  tests.  Equal 
oil  pressure  was  built  up  on  both  sides  of  the  jacks.  In 
order  that  no  load  be  applied  to  the  specimen,  while  the 
pressures  were  being  increased  to  the  required  value,  it 
was  necessary  to  maintain  an  equal  pressure  on  both  sides 
of  the  jacks  at  all  times.  This  was  done  by  keeping 
the  valve  which  interconnects  the  pull  and  push  side  re¬ 
servoirs,  open.  This  valve  was  closed  after  the  required 
pressure  was  attained  on  two  sides  of  the  jack  and  then 
the  pull  side  was  "dumped"  by  opening  the  dump  valve.  The 
latter  was  operated  by  sending  a  predetermined  command 
signal  to  the  servo  valve  which  controls  the  dump  valve. 

In  this  way  the  dynamic  loads  with -a  rise  time  between 
10  to  20  milliseconds  were  applied  to  the  test  specimens. 

A  record  of  all  measurements  such  as  applied  load,  de¬ 
flections,  etc.  with  respect  to  time  was  obtained  on  the 
C.E.C.  recorder  running  at  the  high  speed  of  21.6  inches/sec. 

In  all  cases  a  failure  pulse  slightly  greater 
than  the  estimated  failure  resistance  of  the  specimens 
under  each  type  of  loading, was  applied.  Due  to  limitations 
of  the  loading  machine,  it  was  not  possible  to  apply 
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partial  loads  on  test  specimens  except  for  type  “A" 
specimens  under  compression  mode  loading  (Type  I), 
because  the  failure  loads  were  too  low  in  other  cases. 
Partial  loads  with  a  rise  time  between  10  to  20  milliseconds 
and  a  flat  peak  of  considerably  longer  duration  (about 
1000  milliseconds)  were  applied  to  some  of  these  type  "a'^ 
specimens . 
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REINFORCING  STEEL  -4  NO.  7  WIRES 
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FIGURE  3.1  -  TYPE  A  SPECIMEN 


Figure  3.3  -  Wooden  Form  Work  for  Type  A  Specimens 


Figure  3*4-  "  Device  for  Making  Reinforcement  Cages 
Type  A  Specimen. 
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FIGURE  3.8  -  PUSH  SIDE  RESERVOIR 


HYDRAULIC  SUPPLY  AND  RETURN  LINES 
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FIGURE  3.10- SCHEMATIC  DIAGRAM  OF  LOADING  UNIT 
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Figure  3 *11  “  Experimental  Set  Up  - 
Type  I  Loading. 
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Figure  3*12  -  Experimental  Set  Up  - 
Type  II  Loading. 
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Test  Specimen 
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FIGURE  3.14  -HINGE  SUPPORTS  FOR  TEST  SPECIMEN 
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FIGURE  3.15 -SUPPORTING  DETAIL  "x"FOR  TYPE  A  SPECIMEN 
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FIGURE  3.16 -SUPPORTING  DETAIL  "X  FOR  B  SPECIMEN 
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10  Kips  Load  Cell 
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FIGURE  3.18  LOAD  CELLS  FOR  MEASURING  APPLIED  LOAD 


FIGURE  3.19  -REACTION  LOAD  CELLS 


CHAPTER  1|. 


EXPERIMENTAL  RESULTS 

k-1  INTRODUCTION 

In  this  chapter  the  test  results  of  all  specimens  are 
presented  in  tabular  form  (Tables  4*1  'bo  In  these 

tables  values  of  f  ,  experimental  failure  load,  experl- 

c  ^ 

mental  deflections,  and  rise  time  of  the  failure  load  in 
the  case  of  dynamic  tests  are  given.  Theoretical  failure  - 
load  and  deflections  for  each  specimen,  calculated  on  the- 
basis  of  conventional  theory  presented  in  Section  2-5  of 
Chapter  2,  are  also  given.  It  should  be  noted  that  the 
values  of  the  theoretical  deflections  given  in  each  table 
are  the  elastic  limit  values  based  on  the  failure  load 
obtained  from  the  tests.  Comparisons  between  the  ex¬ 
perimental  results  and  those  obtained  by  the  analysis 
presented  in  Chapter  2  are  made  and  discussed  in  Chapter  5* 
4.2  SUMMARY  OP  TEST  RESULTS  -  TYPE  I  LOADING 
4.2.1  General; 

The  test  results  of  type  ''A"  and  type  "B"  specimens 
are  summarized  in  Tables  4*1  ^nd  4*2,  respectively.  As 
mentioned  in  Section  3*3  of  Chapter  3»  'fche  radial  dis¬ 
tributed  load  was  simulated  by  ten  points  loads  on  the 
arch  specimen.  These  point  loads  are  converted  into 
equivalent  distributed  load  per  inch  and  these  equivalent 
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values  are  presented  in  Tables  I4..I  and  1^.2.  Instead  of 
presenting  the  measured  radial  deflections  at  crown  and  at 
514.°  points  separately,  the  summation  of  these  measured 
deflections  is  presented  in  Tables  14.. 1  and  1^.2  and  compared 
to  the  corresponding  theoretical  values.  This  was  done  in 
order  to  eliminate  unsymmetric  behavior  of  the  test 
specimen,  which  occured  due  to  unavoidable  Irregularities 
in  the  test  specimens  and  in  the  loading  apparatus.  Certain 
specimens  such  as  A-3,  A-9,  A-10  etc.,  or  B-3,  B-I8 
etc.  are  not  included  in  Tables  I4..I  and  }^.2  because  some 
of  these  specimens  were  not  tested  due  to  excessive  honey¬ 
combing  in  the  concrete,  while  for  certain  other  specimens 
reliable  results  were  not  obtained  due  to  malfunctioning 
of  the  loading  device  or  measuring  equipment. 

I4..2.2  Static  Tests; 

The  results  of  the  static  tests  seem  quite 
satisfactory.  The  variation  in  experimental  failure  loads 
is  well  within  +  1S%  of  the  loads  predicted  by  the  conven¬ 
tional  theory  except  for  specimens.  Brl  and  B-2.  In  these 
latter  cases  the  concrete  may  have  had  less  strength  than 
that  Indicated  by  the  test  cylinders.  It  is  also  inter¬ 
esting  to  note  that,  except  for  specimens  A-1  and  B-10, 
the  experimental  failure  loads  are  lower  than  those  pre¬ 
dicted  by  the  conventional  theory.  This  is  probably  due  to 
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secondary  bending  effects  introduced  for  many  reasons.  The 
detailed  discussion  on  the  comparison  between  test  results 
and  the  theoretical  results  is  deferred  until  Chapter  5* 

The  failure  of  the  specimens  was  mainly  by  the 
crushing  of  concrete  at  a  few  points  along  the  arch.  In 
a  few  cases*  the  specimen  failed  by  crushing  at  only  one 
spot,  while  in  other  specimens  the  evidence  of  failure  was 
distributed.  Photographs  of  typical  specimens  after  failure 
are  shown  in  Figs.  I4..I  and  I4..2. 

For  a  few  specimens  such  as  A-1,  A-2,  B-1,  etc., 
reliable  deflection  measurements  were  not  obtained  and 
therefore  these  results  are  not  given  in  Tables  l^.l  and 
Load  deflection  curvesof  typical  specimens  are 
given  in  Fig.  5.I  (Chapter  5)*  Comparison  between 
theoretical  and  experimental  failure  loads  is  shown  in  Fig. 
4.. 3  to  indicate  the  scatter  in  the  test  results. 

4.2.3  Dynamic  Testa: 

In  the  dynamic  teats  of  type  "A"  arches,  specimens 
A-8  and  A-11  were  tested  by  applying  only  one  pulse  corre¬ 
sponding  to  the  failure  pulse.  On  specimens  A-12,  A-13, 
and  A-I7,  a  partial  pulse  was  first  applied  and  measure¬ 
ments  obtained  .  These  specimens  were  then  tested  under  a 
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failure  pulse.  On  specimen  A-19  two  partial  pulses  were 
applied  before  applying  a  failure  pulse.  It  was  not  possible 
to  apply  partial  loads  to  the  type  "B"  specimens  because 
the  failure  loads  are  considerably  lower  and  the  loading 
apparatus  is  not  capable  of  producing  partial  dynamic 
pulses  sufficiently  small. 

Type  "a"  specimens  show  an  average  dynamic  in¬ 
crease  of  about  18^  in  failure  load  over  the  theoretical 
static  failure  loads.  Type  ’'B''  specimens  show  an  average 
increase  of  about  2']%  if  specimens  B-12  and  B-I3  which 
show  exceptionally  large  increases  are  not  considered. 

This  percentage  increase  in  the  failure  loads  seems 
reasonable  on  the  basis  of  the  expected  increase  in  material 
strength.  The  increase  of  about  100^  in  the  failure  loads 
for  specimens.  B-12  and  B-I3  cannot  be  fully  explained 
on  the  basis  of  the  increase  of  material  properties  and 
inertial  effects. 

The  failure  in  all  specimens  was  mainly  by  the 
crushing  of  concrete  simultaneously  at  a  few  points  along 
the  arch.  Photographs  of  typical  specimens  after  failure 
are  shown  in  Pigs.  and  4* 5* 
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The  load-time  curves  and  deflection-time  curves 


of  typical  specimens  are  given  in  Pigs.  5.2,  5.3,  5.4  and  5.5 
of  Chapter  5,  and  the  load-deflection  curve  of  type  *'a" 
specimens  is  shown  in  Fig.  5.6  Chapter  5*  In  this 
curves  the  partial  and  failure  loads  of  each  specimen  are 
non-d Imensionalized  with  respect  to  the  theoretical 
static  failure  loads.  These  non-dimensionalized  values 
are  plotted  against  the  corresponding  experimental 
deflections  values  given  in  Table  4* I*  Tha  experimental 
failure  loads  are  plotted  vs.  theoretical  static  loads 
in  Fig.  4*3  to  Indicate  the  scatter  in  test  results  and 
the  general  difference  between  static  and  dynamic  results. 

4.3  SUMMARY  OF  TEST  RESULTS  -  TYPE  II  LOADING 
4«3*1  General: 

The  test  results  for  type  ''A”  and  type  ''B" 
specimens  are  summarized  in  Tables  4*3  and  4*4,  respectively. 
4.3*2  Static  Tests: 

The  results  of  the  static  tests  seam  quite  reason¬ 
able.  The  failure  load  of  B-19  is  within  2.%  of  the  failure 
load  predicted  by  the  conventional  theory.  While  for 
specimens  B-21,  A-21,  and  A-22,  the  experimental  values 
are  higher  by  10  to  15^*  There  are  several  factors  which 
could  have  caused  the  higher  experimental  values.  One 
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reason  could  be  that  the  supports  may  not  behave  like 
perfect  hinges  and  may  introduce  some  restraint,  which  would 
tend  to  increase  the  experimental  failure  load.  Secondly, 
the  yield  stress  f  in  the  steel  is  taken  at  ij.0  ksi, 
while  predicting  the  conventional  theoretical  loads.  The 
stress-strain  curve  of  the  steel  (See  Pigs.  3*5  and  3*6) 
is  not  exactly  bi-linear  but  follows  a  flat  curve  after 
[[.0  ksi.  It  is  obvious  that  the  strain  in  the  steel  at 
failure  in  both  ''A"  and  ''B"  specimens  is  much  greater 
than  the  yield  strain  at  4.0  ksi,  which  suggests  that 
at  failure  the  stress  could  be  higher  than  4^  ksi. 

Therefore,  the  experimental  failure  loads  could  be  higher 
than  that  predicted  by  conventional  theory  using 
f  =40  ksi. 

y 

The  failure  of  the  specimens  occured  by  the 
formation  of  three  hinges,  one  at  the  crown  and  the 
other  two  in  the  vicinity  of  the  quarter  points.  The 
locations  of  these  latter  hinges  almost  coincided  with  the 
theoretical  prediction  of  hinge  formation  at  points  38° 
from  the  supports.  Photographs  of  typical  specimens 
after  failure  are  shown  in  Figs.  4*6  and  4*7« 
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The  failure  of  the  specimens  seemed  quite  ductile. 
The  maximum  radial  deflection  at  the  crovm  was  of  the 
order  of  0.1|5"  for  type  *'a"  and  0.85“  for  type  specimens. 

The  load-deflection  curves  of  typical  specimens 
are  given  in  Fig.  5.7  and  5.8  Chapter  5.  To  represent  the 
scatter  in  test  results  the  experimental  failure  loads 
are  plotted  against  conventional  theoretical  loads  in 
Pig.  I4..8. 

4*3«3  Dynamic  Tests: 

All  specimens  were  tested  for  failure  only. 

Because  the  failure  loads  were  extremely  small,  no 
attempt  was  made  to  apply  partial  loads. 

Type  "a'"  specimens  show  an  average  dynamic 
increase  of  about  55l^  'the  failure  loads  over  con¬ 
ventional  theoretical  static  failure  loads.  For 
Type  "B"  specimens  the  increase  is  of  the  order  of  70^. 

Part  of  this  rather  large  dynamic  increase  can  be 
attributed  to  the  larger  than  expected  static  strength 
as  discussed  in  the  preceeding  section.  The  remaining 
increase  is  due  to  the  increase  in  material  properties 
due  to  very  rapid  strain  rates  and  to  inertial  effects. 

On  this  basis,  the  test  results  seem  quite  reasonable. 
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The  failure  of  the  specimens  was  similar 
to  that  observed  in  static  tests  as  described  in 
Section  4-3 -2,  Photographs  of  typical  specimens  after 
failure  are  shown  in  Figs.  4.9  and  4.10. 

The  failure  of  the  specimens  seemed  quite 
ductile.  The  maximum  radial  deflection  at  the  crown 
was  about  0.3"  for  type  A  and  about  0,9  for  type  B 
specimens . 

The  load-time  and  deflection-time  curves  of 
typical  specimens  are  shown  in  Pigs.  3.9  and  3.10 
Chapter  3.  The  scatter  in  test  results  is  indicated 
by  Pig.  4.8. 

4.4  SUMMARY  OP  TEST  RESULTS  -  TYPE  III  LOADING 
4.4.1  General: 

The  test  results  for  type  A  and  type  B  specimens 
are  summarized  in  Tables  4.3  and  4.6,  respectively.  The 
results  for  specimen  B-28  are  not  given  in  Table  4-6, 
because  no  reliable  results  were  obtained. 

The  deflection  data  in  Tables  4.3  and  4.6,  clearly 
indicate  that  the  strain  in  the  steel  wires  at  failure 
must  have  been  considerably  higher  than  the  strain 
corresponding  to  a  stress  of  4^  ksi,  at  which  yielding 
begins.  Therefore,  as  explained  in  Section  4.3.2,  it 
seems  reasonable  to  assume  a  higher  value  of  fy  while 
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calculating  the  failure  load  based  on  conventional  theory. 

A  value  of  I4.6  ksl  Is  used  In  the  calculation  of  these 
conventional  theoretical  failure  loads. 

4..[|..2  Static  Tests: 

The  results  of  the  static  tests  of  type  "a" 
specimens  seem  quite  reasonable.  Experimental  failure 
loads  on  an  average  are  about  9^  higher  than  the 
conventional  theoretical  failure  loads,  except  specimen 
A-33j  for  which  the  failure  load  Is  lower  by  about  95?. 

The  latter  specimen  was  actually  loaded:  twice.  The  first 
run  was  terminated  after  about  60^  of  the  failure  load 
had  been  applied,  due  to  difficulties  In  the  loading 
device.  It  was  then  loaded  to  failure.  The  preloading 
could  have  damaged  the  specimen  and  caused  failure  under 
a  lower  load.  The  increase  in  the  failure  load  of  all  other 
specimens  over  the  conventional  theoretical  load  seems 
partly  due  to  the  restraint  at  the  support  as  explained 
in  Section  I4..3.2  and  partly  due  to  the  effect  of  the 
weight  of  the  loading  jacks  on  the  test  specimen  as 
explained  below. 

It  Is  clear  from  Pig.  3*13  that  about  half  the 
weight  of  the  jack  la  applied  to  the  arch  at  the  loading 
points.  Theoretical  considerations  Indicate  that  if  the 
top  jack  was  "pulling"  the  arch  and  the  bottom  jack  "pushing" 
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(see  Pig.  3.13) >  the  failiwe  load  increases  by  the  same 
amount  as  the  applied  weight  of  the  jack.  On  the  other 
hand,  with  the  top  jack  "pushing"  and  the  bottom  jack 
"pulling",  the  failure  load  decreases  by  the  same  amount. 

This  effect  is  more  severe  in  the  case  of  type 

B  specimens  because  the  failure  loads  are  only  six  to 
seven  times  the  applied  weight  of  the  jack  (which  is  about 

10  lbs)  as  compared  to  I4.O  times  the  applied  weight  in  the 

case  of  type  A  specimens.  In  order  to  eliminate  this 

effect,  specimens  B-26  and  B-27  were  tested  with  the  top  jack 

"pulling"  and  the  bottom  jack  "pushing",  while  the  reverse 

was  done  for  B-29  and  B-30.  Table  4.6  indicates  that  the 

average  of  the  failure  loads  of  B-26  and  B-27  is  91  lbs, 

and  of  B-29  and  B-30  is  52  lbs.  Therefore,  an  average  of 

these  two  values  or  about  71  lbs.  should  be  considered  as 

the  experimental  failure  load.  The  average  conventional 

theoretical  load  for  these  specimens  is  about  60  lbs.,  and 

the  eleven-pound  increase  in  the  experimental  load  can 

easily  be  due  to  the  various  factors  mentioned  previously, 

or  small  errors  in  measurement.  The  test  results  are 

therefore  considered  quite  reasonable. 

The  failure  of  the  specimens  occurred  by  the  form¬ 
ation  of  hinges  very  nearly  at  the  qT;iarter  points.  Photo- 


-116- 


graphs  of  typical  specimens  after  failwe  are  shown  in 
Pigs.  I4..II  and  I4..12.  The  failures  were  quite  ductile  with 
maximum  radial  deflection  at  the  quarter  points  of  the  order 
of  1.5"  for  type  A  specimens  and  about  2"  for  type  B 
specimens . 

The  load  deflection  curves  of  typical  specimens, 
are  given  in  Pigs,  5.11,  5.12  and  5.13,  Chapter  5.  The 
experimental  failure  loads  are  plotted  against  conventional 
theoretical  loads  in  Pig.  4*13  to  indicate  the  scatter  in 
experimental  results. 

4.4.3  Dynamic  Tests: 

All  specimens  were  tested  for  failure  only. 

No  attempt  was  made  to  apply  partial  loads  for  the  reasons 
already  mentioned. 

Type  A  specimens  show  an  average  increase  of 
about  Q%  in  failure  loads  over  the  conventional  theoretical 
static  failure  loads.  Out  of  this  about  10^  is  probably 
due  to  partial  support  restraint  and  the  effect  of  the 
weight  of  jacks  on  the  failure  loads,  as  already  explained 
in  previous  sections.  The  remaining  increase  is  attributed 
to  the  increase  in  the  material  properties  due  to  rapid 
strain  rates  and  Inertial  effects.  This  is  discussed  in 
more  detail  in  the  next  chapter. 
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Dynamic  tests  on  type  B  specimens  were  performed 
with  the  top  jack  "pulling"  and  the  bottom  jack  "pushing". 

No  tests  were  performed  in  the  reverse  condition  as  was 
done  in  static  tests.  Therefore,  the  results  of  dynamic 
tests  of  these  specimens  should  be  compared  with  the  corres¬ 
ponding  results  of  the  static  tests,  i.e,,  the  results  for 
specimen  B-26  and  B-27.  This  type  of  comparison  shows  an 
average  Increase  of  l80^  in  dynamic  failure  loads.  This 
increase  seems  to  be  due  to  the  increase  in  the  material 
properties  under  very  rapid  strain  rates  and  inertial 
effects.  This  aspect  is  discussed  further  in  the  next 
chapter. 

The  failure  of  the  specimen  occurred  by  the 
formation  of  hinges  very  nearly  at  the  quarter  points,  as 
in  the  static  tests.  Photographs  of  typical  specimens 
after  failure  are  shown  in  Pigs,  4.11+  and  1^,1$.  The  failure 
of  the  specimen  was  quite  ductile.  The  maximum  radial 
deflection  at  quarter  points  is  about  1.25"  for  type  A  and 
1.7"  for  type  B  specimens. 

The  lead-time  and  deflection-time  curves  for 
typical  specimens  are  shown  in  Pig.  5.11+  and  5.15  in 
Chapter  5.  To  Indicate  the  scatter  of  experimental  results, 
experimental  failure  loads  are  plotted  against  theoretical 
static  failure  loads  in  Pig.  1+.13. 
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TABLE  4.1 


SUMMARY  OF  TYPE  "A"  SPECIMENS  -  TYPE  I,  LOADING 
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These  are  loads  representing  flat  peak  of  the  partial  pulse. 

This  is  the  summation  of  radial  deflection  at  crown  and  at  one  5I|.®  point 

These  are  the  summations  of  max.  deflection  at  crown  and  at  54°  points,  at 
corresponding  to  the  flat  peak  at  partial  pulse. 

-  Experimental  failure  load. 

THgT  _  Theoretical  static  failure  load. 
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-  Summation  of  maximum  experimental  radial  deflections  at  crown  and  at  two  54°  points. 


TH 


Summation  of  theoretical  elastic  radial  deflection  at  crown  and  at  two  54°  points, 
at  experimental  failure  load. 


t 
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-  Rise  time  of  the  load. 
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OF  TXPE  "B"  SPECIMENS  -  TYPE  I  LOADING 
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"Notations  use 


OP  TYPE  "A”  SPECIMENS  -  TYPE  II  LOADING 
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Rise  time  of  the  load 
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SUMMARY  OF  TYPE  "B"  SPECIMEN  -  TYPE  III-  LOADING 


1+.5'  NATURAL  PERIOD  OF  ARCH  SPECIMENS 

The  natural  periods  of  the  arch  for  the  first  bending 
or  antisymmetric  mode  and  the  second  bending  mode  (i.e., 
the  mode  corresponding  to  the  configuration  of  the  arch, 
under  a  point  load  at  the  crown)  were  experimentally 
determined  for  one  specimen  of  each  type.  These  values 
are  given  in  Table  I)..?.  The  natural  period  for  the 
compression  mode"  was  not  determined  experimentally. 


TABLE  4.7 
NATURAL  PERIODS 


1 

1  Mode 

1 

Natural  Period 

■ 

(Milliseconds)  | 

j 

» 

1 

Type  A 

Type  B 

! 

( 

Specimen 

Specimen 

First  bending  mode 

25 

59 

Second  bending  mode 

6 

12 

"The  natural  period  in  the  compression  mode  is  approximately 
given  by  r/l800  (see  Ref.  No.  13  in  bibliography,  p.  5B-15)* 
Therefore  the  natural  period  for  both  types  in  the  compression  mode 
is  about  0.83  milliseconds.  In  the  above  formula  "r"  is  the 
radius  in  feet. 
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Figure  [4..I  -  Appearance  of  Specimen  A-I4.  after  Test 
(Static  Test  -  Type  I  Loading). 
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FIGURE  4.3-  p,,  VS  p,,  -TYPE  I  LOADING 
f'UEXP  “THst 
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Figure  4-. 6  -  Appearance  of  Specimen  A-21  after  Test  - 
(Static  Test  -  Type  II  Loading). 
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Figure  I4..7  -  Appearance  of  Specimen  B-21  after  Test  - 
(static  Test  -  Type  II  Loading). 
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FIGURE  4.8-  Pugj^p  VS.  TYPE  H  LOADING 
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Figure  I4..II  -  Appearance  of  Specitnen  A-30  after  Test 
(Static  Test  -  Type  III  Loading). 
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Figure  Ix.lZ  -  Appearance  of  Specimen  B-26  after  Test  - 
(Static  Test  -  Type  III  Loading) 
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Pi.gure  Ij-ilLj.  -  Appearance  of  Specirnen  A-  35  after  Teat  - 
(Dynamic  Test  -  Type  III  Loading). 
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CHAPTER  5 


COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
RESULTS 

5.1  INTRODUCTION 

For  the  purpose  of  comparison,  certain  experimental 
specimens  of  both  type  A  and  type  B  are  analysed  theoretically. 
The  nonlinear  theory  developed  In  2.1  through  2.I4.  Is  used 
to  obtain  the  static  and  dynamic  response,  ultimate  loads  and 
/I's  for  type  I  and  type  III  loading.  The  average  dynamic 
strain  rate  6(3^  necessary  to  obtain  the  dynamic  properties  of 
concrete  (equation  2.2b)  is  calculated  as  e^c  =  eu/tr*  These 
ultimate  loads  and  are  called  analytical  results.  The 

elastic  limit  deflections  needed  to  calculate  the  ^*s  are 
obtained  by  using  the  formulae  presented  in  ^2.$.  Further, 
a  set  of  static  ultimate  loads  for  all  three  types  of  loading 
is  also  obtained  from  the  conventional  theory  presented  in  ^2.5* 
The  ultimate  loads  calculated  in  this  way  are  referred  to  as 
"theoretical  ultimate  loads". 

5.2  THEORETICAL  AND  EXPERIMENTAL  RESULTS  AND  THEIR 

COMPARISON. 

5.2.1  Type  I  Loading; 

A  summary  of  the  experimental  and  analytical 
values  of  the  ultimate  loads  and  yu’s  for  certain  specimens 
of  type  A  and  type  B  is  given  in  tables  S.la  and  ^.Ih 


-141- 


respectively.  Also  listed  in  the  tables  are  the  theoretical 
static  ultimate  loads. 

(a)  Static  Behavior: 

Tables  5«la  and  ^.Ih  clearly  show  that 
the  analytical  values  of  the  ultimate  loads  agree  within 
ten  percent  with  the  experimental  values  for  specimens  A-2, 
A-l+j  A-5,  B-11,  B-15,  and  B-16.  However,  the  experimental 
ultimate  loads  for  specimens  B-1  and  B-2  are  much  lower  and 
those  for  A-1  and  B-10  are  higher  than  the  analytical  values. 
This  may  be  due  to  a  possible  difference  between  the  strength 
of  concrete  in  the  test  specimens  and  that  Indicated  by  the 
test  cylinders. 

The  comparison  between  the  theoretical  and  experi¬ 
mental  ultimate  loads  has  already  been  made  in  I4..2.2.  How¬ 
ever,  attention  la  drawn  here  to  the  fact  that  the  experiment¬ 
al  failure  loads  are  in  practically  all  cases  lower  than  the 
theoretical  failure  loads.  This  seems  to  be  due  to  the 
presence  of  secondary  bending  effects  which  are  not  considered 
in  the  theoretical  calculations.  The  secondary  bending 
effects  occur  for  such  reasons  as:  (1)  the  shape  of  the  arch 
does  not  confirm  to  the  funicular  polygon  of  the  compression 
mode  loading;  (2)  in  the  tests,  the  uniformly  distributed 
load  is  simulated  by  ten  equally  spaced  point  loads; { 3)These 
ten  point  loads,  applied  on  the  test  specimen  may  not  be 
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exactly  equal,  due  to  the  limitations  of  the  loading 
apparatus;  ([].)  the  arch  cross-section  may  not  be  exactly 
uniform  along  the  length. 

A  comparison  between  the  analytical  and  theoretical 
values  shows  that  both  agree  within  ten  percent  for  most  of 
the  specimens.  The  analytical  values  are,  however,  consist¬ 
ently  lower  in  all  cases;  the  reason  being  that  the  secondary 
bending  effects  are  not  considered  in  calculating  theoretical 
loads . 

Table  $.1  (a)  and  5*1  (b)  also  Indicate  that  except 
for  specimen  A-5»  the  experimental  -  values  are  lower  than 
the  analytical  m  -  values. 

Fig.  5*1  shows  the  analytical  and  experimental  load- 
deflection  curves  for  specimen  A-5»  selected  for  illustration, 
to  be  in  fair  agreement. 

(b)  Dynamic  Behavior: 

It  is  seen  from  tables  S.la  and  that  the 

analytical  values  of  dynamic  ultimate  load  agree  within 
about  fifteen  percent  with  the  experimental  values  for  all 
specimens  except  specimens  B-12  and  B-20.  The  high  experi¬ 
mental  ultimate  loads  for  these  specimens  may  be  in  part  due 
to  an  Increase  in  material  strength  under  rapid  strain  rate, 
higher  than  that  considered  in  the  analytical  approach. 

The  analytical  and  experimental  ju-values  lie 
within  a  range  of  1.^  to  2.5- 
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Figures  5*2  5»3  and  5*5  show  analytical  and 
experimental  load-time  and  deflection-time  curves  for 
specimens  A-11,  A-l?  and  B-II4..  These  curves  represent  a 
failure  load  pulse.  Specimen  A-l?  was  also  tested  under 
partial  load  pulse  and  the  load-time  and  deflection-time 
curves  for  this  partial  load  pulse  are  shown  in  figure  5*4* 
All  the  curves  show  a  fair  agreement  between  the  analytical 
and  experimental  results.  Figure  5*6  shows  the  analytical 
and  experimental  dynamic  load-deflection  curve  for  type  A 
specimens.  The  points  on  the  experimental  curve  are  obtained 
from  a  number  of  tests  by  nondimensionalizing  the  dynamic 
load  with  respect  to  the  corresponding  theoretical  static 
ultimate  load  of  each  test.  The  analytical  curve  pertains 
to  specimen  A-l?. 

(c)  Comparison  Between  Dynamic  and  Static 

Behavior; 

Table  5*3  shows  the  experimental  and 
analytical  values  of  the  dynamic  Increase  factors  (DIF)  for 
both  type  A  and  type  B  specimens.  The  dynamic  Increase 
factor  is  obtained  as  a  ratio  of  the  average  nondlmensional 
dynamic  load  to  the  nondlmensional  static  load;  the  average 
being  obtained  from  the  values  given  in  tables  and 

The  values  of  experimental  ultimate  loads  of  specimens  B-12 
and  B-I3  are  not  considered  in  obtaining  the  average  non- 
dimensional  dynamic  load  because  these  values  seem 
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unreasonably  high.  The  experimental  and  analytical  dynamic 
Increase  factors  vary  between  1.2  and  l.I|..  This  increase 
is  mainly  due  to  the  Increase  in  the  properties  of  materials 

at  very  rapid  strain  rates. 

5.2.2  Type  II  Loading 

No  comparison  between  the  experimental  and 
analytical  results  under  this  type  of  loading  is  possible 
because  the  non-linear  theory  cannot  be  used  to  predict  the 
ultimate  loads  for  the  reasons  explained  in  ^  2.1.3*  A 
comparison  between  the  experimental  and  theoretical  failure 
loads  is  already  made  and  reasons  for  the  differences  in  the 
values,  if  any,  are  explained  in  ^  14..3.2.  Therefore,  in 
what  follows,  only  a  comparison  between  the  experimental 
static  and  dynamic  behavior  is  made. 

Ih  table  5*4>  'cne  values  of  static  and 
dynamic  ultimate  loads  for  the  type  A  and  type  B  arches  are 
given.  These  are  the  average  values  of  those  given  in 
tables  [|..3  and  I4..4  for  type  A  and  type  B  specimens  respec¬ 
tively.  It  is  reasonable  to  take  such  an  average,  even 
though  the  value  of  f^  varies,  for  each  specimen  in  the 
testa,  because  the  Influence  of  f^  on  the  ultimate  loads 
of  such  underreinforced  elements,  where  the  bending 
moment  is  predominant  compared  to  the  axial  thrust,  is 
insignificant.  Table  5*4  shows  that  the  dynamic  increase 


is  of  the  order  of  36/C  for  type  A  arch  and  about  37^  for 
type  B  arch.  This  increase  seems  to  be  partly  due  to  the 
increase  in  the  material  properties  under  very  rapid  strain 
rates  and  partly  due  to  inertial  effects. 

The  load-deflection  curves  for  specimens  A-21 
and  B-21  (static  tests)  are  shovm  in  Figures  5*7  and  5*8» 
while  Figures  5*9  and  $-10  show  the  load-time  and  deflection 
time  curves  for  specimens  A-2i4.  and  B-2I4.  (dynamic  tests). 

3.2.3  Type  HI  Loading 

The  experimental  values  of  the  static  and 
dynamic  ultimate  loads  and  ^'s  for  type  A  and  type  B  arches 
are  given  in  table  ^.2.  The  values  of  the  theoretical 
static  ultimate  loads  are  shown  in  table  ^,2.  These  values 
are  the  averages  of*  those  presented  in  tables  and 
The  analytical  values  of  the  static  ultimate  loads  and  /I's 
given  in  table  5*2  are  those  obtained  for  specimens  A-3I1- 
and  B-26,  while  the  corresponding  values  for  dynamic  tests 
are  for  specimens  A-37  and  B-32.  The  analytical  values  are 
calculated  by  taking  fy  =  I4.6  ksi,  for  the  reasons  explained 
in  article  I4..I4..I,  A  comparison  between  these  analytical 
values  and  the  average  experimental  values  is  reasonable, 
for  reasons  already  mentioned  in  the  previous  section. 

(a)  Static  Behavior: 

The  experimental  ultimate  loads  agree 
with  the  analytical  and  theoretical  values  within  about 


-146- 


ten  percent.  Table  5*2  also  indicates  an  excellent 
agreement  between  theoretical  and  analytical  values  of 
ultimate  loads.  The  analytical  values  of  ^'s  are  much 
smaller  thsin  the  experimental  Ja^  s  because  of  the  convergence 
difficulties  in  the  analytical  approach  as  explained  in 
f  2.3.7.  Figures  ^•ll,  5*12  and  5*13  show  the  analytical 
and  the  experimental  load  deflection  curve  for  specimens 
A-28,  A-3I4-  and  B-26  respectively. 

(b)  Dynamic  Behavior: 

The  analytical  value  of  the  ultimate  load 
of  the  type  A  specimen  is  17%  lower  than  the  experimental 
value  while  for  the  type  B  specimen,  this  difference  is  as 
high  as  4-0^.  The  analytical  ^‘s  are  lower  than  the 
experimental  /a's.  The  lower  analytical  u-value  can  be 
explained  on  the  basis  that  in  the  analytical  calculations, 
the  arch  is  considered  to  have  failed  when  the  compressive 
strain  in  extreme  fibres  at  any  one  point  in  the  arch 
exceeds  e^  =  O.OO38  inches/  in.:  while  in  tests,  the 
extreme  fibre  strains  at  failure  could  be  much  larger  than 
e^.  In  figures  5.14  5* 3-5  are  shown  the  analytical  and 

experimental  load-time  and  deflection-time  curves  for 
specimens  A-37  and  B-32  respectively.  These  curves  indicate 
a  fair  agreement  between  the  analytical  and  experimental 
results. 
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(c)  Gompariaon  of  Dynamic  and  Static  Behavior: 

In  table  a  comparison  between  the 

experimental  and  analytical  dynamic  increase  factor  (DIF) 
is  sought.  The  value  of  the  experimental  static  ultimate 
load  for  type  B  arches,  given  in  this  table  is  an  average 
of  those  of  specimens  B-26  and  B-27j  while  the  specimens  B-29 
and  B-30  are  not  considered  in  thid  average  value,  in  order 
to  have  a  consistent  comparison  between  static  and  dynamic 
behavior  as  explained  in  detail  in  paragraph  3  of  ^  4-*4-*3* 

The  experimental  values  of  the  dynamic  increase 
are  of  the  order  of  65^  and  179/^  respectively  for  the  type  A 
and  type  B  arches.  The  corresponding  analytical  values  are 
and  Thus  the  agreement  between  the  experimental 

and  analytical  values  is  quite  good.  This  dynamic  Increase 
seems  to  be  due  to  the  increase  in  the  material  properties 
under  very  rapid  strain  rates  and  also  due  to  Inertial 
effects.  The  comparatively  larger  dynamic  Increase  in  the 
case  of  type  B  arches  is  possibly  due  to  the  presence  of 
substantial  inertial  effects, as  compared  to  type  A  arches. 

The  experimental  as  well  as  analytical  _>u-value3 
for  the  dynamic  case  are  smaller  than  those  for  the  static 
case.  (See  table  5*2).  This  is,  in  part,  due  to  the  fact 
that  the  elastic-limit  deflections  are  calculated  at  the 
ultimate  loads  which  are  considerably  higher  in  the  dynamic 


case.  Therefore  the  elastic-limit  deflections  are 
considerably  higher  and  the  /i-values  are  lower. 
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TABLE  g.la" 

COMPARISON'  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  TYPE  I  LOADING 
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TABLE  5.3 


COMPARISON  OF  STATIC  AND  DYNAMIC  BEHAVIOR 
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FIGURE  5.1-SPECIMEN 
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FIGURE  5.2-SPECIMEN  A-ll  {FAILURE  PULSE) 
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ANALYTICAL  RESULTS 
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FIGURE  5.15-SPECIMEN  B-32 


CHAPTER  ^ 


CONCLUSIONS 

The  following  conclusions  can  be  drawn  on  the  basis 
of  the  analytical  and  experimental  investigations  as  well 
as  the  approximate  theoretical  analysis. 

1)  The  static  ultimate  loads  obtained  experimentally 
are  in  good  agreement  with  both  the  conventional 
approximate  theory  and  the  non-linear  theory  developed 
herein. 

2)  In  the  compression  mode  loading,  the  dynamic 
increase  in  ultimate  load  is  mainly  due  to  the 
increase  in  properties  of  materials  under  very 
rapid  strain  rates.  This  increase  is  about  30 
to  1%, 

3)  In  the  eases  of  a  concentrated  load  at  the 
crown  and  antisymmetric  concentrated  loads  at 
quarter  points,  the  natural  periods  of  the  arches 

in  the  flexural  mode  (and  hence  the  inertial  effects) 
have  a  significant  Influence  on  the  dynamic  Increase. 
This  Increase  is  much  higher  than  that  in  the  case  of 
compression  mode  loading. 

4.)  In  the  compression  mode  loading  the  analytical 
and  experimental  values  of  p.  are  of  the  order  of 

1.5  to  2.5. 
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5)  For  a  concentrated  load  at  the  crown  the 
experimental  value  of  jx  la  about  9  fbr  the 
static  case  and  about  7  for  the  dynamic  case. 

6)  The  experimental  value  of  jx  for  antisymmetric 
concentrated  loads  at  quarter  points  is  of  the 
order  of  1$  for  the  static  case  and  varies  between 
3  and  7  for  the  dynamic  case. 

7)  The  analytical  approach  is  adequate  to  compute 
the  ultimate  loads  in  moat  of  the  cases  considered; 
however,  it  gives  quite  conservative  results  for 
the  dynamic  antisymmetric  quarter  point  loading. 

8)  Experimental  and  analytical  results  indicate 
that  the  approximate  static  theory  is  quite  adequate 
to  predict  the  static  ultimate  loads  for  under- 
reinforced  sections.  The  d3mamlo  ultimate  loads 
for  compression  mode  loading  can  also  be  predicted 
by  the  approximate  theory  provided  that  an  approp¬ 
riate  dynamic  Increase  factor  (based  on  the  increase 
in  material  properties)  is  used. 

9)  Prom  the  comparison  between  the  analytical 
results  and  those  obtained  from  the  approximate 
analysis  based  on  a  single  degree  freedom  system 
(  ^  2.7)  it  seems  that  the  latter  may  give 
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unconservative  dynamic  failure  loads,  if  the  duration 
of  loading  is  less  than  about  half  the  natural  period. 
However,  clear  reasons  to  explain  this  fact  are  not 
evident  and  further  investigation  is  needed. 
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APPENDIX  I 


VARIOUS  CONSIDERATIONS.  FOR  FORCE-STRAIN  RELATIONS 
BASED  ON  LINEAR  CONCRETE  STRESS-STRAIN  CURVE 


As  explained  in  §  2.3.I(.,  in  order  to  solve  the  non¬ 
linear  force-strain  relations  to  obtain  strains  ei  and 
e^,  it  is  necessary  to  have  certain  approximate  values 
of  these  strains  to  start  the  ‘Newton-Raphson  iteration* 
(Appendix  II).  These  approximate  values  of  strains  at 
the  bottom  and  the  top  of  the  section  denoted  here  by  eei 
and  ee4  are  calculated  by  assuming  that  the  stress-strain 
curve  for  concrete  is  linear  with  a  modulus  of  elasticity 
E^,  and  that  the  tensile  stresses  in  concrete  are  neg¬ 
ligible.  The  value  ^of  E^  is  assumed  to  be  less  than 
the  initial  tangent  modulus,  so  that  the  linear  stress- 
strain  curve  closely  approximates,  the  non-llnear  curve 
even  at  high  stress  levels.  The  following  considerations 
are  necessary  to  calculate  the  values  of  eei  and  ee^. 

At  any  section  1,  Nj^  and  are  known.  Prom  these 
are  calculated. 


II  II  B 

N  =  N  .  f„bt  and  M  =  M  f  „bt 

C  V 


(Al.l) 


Also  are  calculated  the  gross  Area  A  and  the  gross  moment 

S 

of  inertia  I  of  the  section, 

g 
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{A1.2) 


and 


A  =  bt  +  p,bt(n-l) 

Ig  =-^  +  p^bt(n-l) 


The  various  possibilities  of  the  strain  distributions 
on  the  section  are, 

{i)  N  =  0  and  M  is  positive  or  negative; 


Figure  Al.l  Strain  Distribution  -  Case  (1) 


(ii)  N  is  positive  and  M  is  positive  or  negative; 

There  exist  two  possibilities  in  this  case. 
2  X 

If  )  the  entire  section  is  under 

N  tA 

^  g 

compression  (Fig.  A1.2);  otherwise  the  section  is  partly 
in  compression  and  partly  in  tension  (Pig.  A1.3). 


Figure  A1.2  Strain  Distribution  -  Case  (11) 
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Figure  AI.3  Strain  Distribution  -  Case  (ii  ) 


Figures  AI.I4.  Strain  Distribution  -  Case  (iii) 


(iii)  N  is  negative  and  M  is  positive  or  negative: 

Again  there  exist  two  possibilities.  The  first 
possibility  is  that  the  section  is  partly  in  compression 
and  partly  in  tension  (Pig.  Al.tj),  while  the  second 
possibility  is  that  the  section  is  entirely  under  tension 
(Fig.  A1.5)»  The  criterion  used  for  distinguishing  these 
possibilities  from  each  other  is: 

If  ^  0,  Pig.  AI.5  is  applicable,  and 

if  otherwise,  Fig.  Al.  ij.  is  to  be  used.  This  criterion 
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indicates  that 


In!  M  / 

is  approximate  since  0 

the  strains  in  the  concrete  cover  over  the  top  steel  ( if 
M  is  positive) or  in  the  concrete  cover  over  the  bottom 
steel  ( if  M  is  negstive)  are  compressive.  This  signifies 
that  the  section  is  not  completely  in  tension  until  |N| 
is  somewhat  greater  than  that  required  by  this  criterion. 
However,  since  the  depths  of  concrete  covers  are  usually 
small  as  compared  to  the  total  depth  of  the  section,  this 
approximate  criterion  is  accepted. 


-ee 

I 


^ ^ 


Figure  Al« 5  Strain  Distribution  -  Case  (lii) 

The  equations  necessary  to  calculate  the  strains 
ee  and  ee  in  all  the  above  cases  are  readily  obtained 

1  4 

by  considering  the  equllibriuni  of  the  internal  and 
external  forces  acting  on  the  section.  These  equations 
^nd  the  criteria  used  to  classify  all  the  cases  are 
incorporated  In  the  computer  program,  (Appendix  HI). 
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APPENDIX  II 


NONLINEAR  SIMULTANEOUS  ALGEBRAIC  EQUATIONS 


A2.1  'NEWTON-RAPHSON'  METHOD 
( 12) 

This  method  affords  an  effective  iterative  procedure 

to  obtain  the  solution  of  two  or  more  simultaneous  non¬ 
linear  equations.  For  example,  the  two  simultaneous 
equations  of  the  form, 

f(x,y)  =  0  ,  g{x,y)  =  0  (A2.1) 


having  (cC,  /5)  as  one  of  its  real  solutions  can  be  solved 
to  obtain  oC  and  as  follows: 

The  two  functions  can  be  expanded  in  a  Taylor  series 


0  =  f(«c,/5)  = 

“  +  + . 


where  and  y^^  are  the  values  of  the  iteration,  If 

OC  and  ^  in  the  terms  on  the  right-hand  side  of  the 


-177- 


The  iteration  la  started  at  a  value  (x  ,  y  )  which 

O  0 

is  sufficiently  near  )  and  equations  (A2.3)  and 

(A2.1;)  are  used  to  obtain  successively  better  values  of 

the  roots  of  the  equations  (A2.1).  The  iteration  is 

continued  until  the  corrections  become  smaller  than  some 

assigned  tolerance  limit.  The  advantage  of  this  method 

over  others  is  that  this  is  a  *  second-order’  processj  that  is, 

,  .th 

when  the  iteration  converges,  the  errors  in  the  (k+1) 
iterate  tend  to  be  linear  combinations  of  the  'squares' 
of  the  errors  in  the  iterate. 
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A2.2  APPLICATION  OP  THE  METHOD  TO  THE  PRESENT  PROBLEM 
As  is  explained  in  §  2.3*4»  this  method  is  used  to 
solve  the  governing  equations  (2.26)  and  (2.27)  or  (2.28) 
and  (2.29).  The  forma  of  equations  used  are  as  follows; 

Using  the  notation 

X  =  ei  ,  y  =  64  , 

k  =  ^t  ,  (1-r)  and  ki  =  ^s  ^t 
— n —  n  - n — 

2f  2f 

c  c 

and  substituting 

eg  =  e 1  + 

and  03  =  61  + 


(64  ■■  ®  1) 
(64  “  ®  1) 


(t-d) 

t 

d 

t 


equations  (2.2'6}  and  (2*27)  become, 
f(x,y)  =  Wx*y)  X  iJL-p.  . 


-  N  =  0 


3® 


g(x,y)  -  2  - 


_L  (^-x)^  + 

a  vyx/T-  , 

^  6e^  12e  ^ 


>  (A2.5) 


-  M(y-x)  =  O' 

when  62  e^_  and  63  ^  e^. 
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Similarly  equations  (2*28)  and  (2.29)  become 


f(x,y)  =  ki  (y  -X  )  -  --  xy 


Ic  1  k  1  ^  1 


il  ..  -j1 

I  1 2 

e  3e„ 

c  c 


g(x,  y)  - 


2  t^ 


(y-x)' 


-  N  (y-x)  =  0 


A-  x(y-x) 
2  t  n 


HA2.8) 


3 

+  L 


6e, 


12e, 


+ 

2e. 


y°x 

6e'' 


M  (y-x)  =  0 


when 


e.  ^  ey  and  e,  ^ 


J 

ki  ki  ^  2^ 


f(x,y)  =  -  ke^  (y-x)  +  (ki  -  ^  +(ki--  — ) 

+  Z._  _  iL_2  -  N  (y-x)  =  0 

e„  3e^ 


g(x,y)  =  (1-i)  x(y-x)^.^  l-'-f  (1-  J)(y-x) 


ki  ' 


f  e^.  (y-x)^  .  A  - 

2  ^  y  6e  12e„ 


+  2_4_  _  M  (y-x)^  =  0 

6e„ 


s 

iUL 

I 

2e„ 


>{  A2.9) 


when  1 68 1  ^  |e^l  and 
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f  {x,y) 


g(x,y) 


when 


kix  (y-x)  +  ki  — -  (y-x)®  +  kie  (y-x) 

^  y 


x(y-x)  -  ~  ^  (y-x)=  + 

e 


n 


3e, 


*  2 


N  (y-x)  =  0 


d*  , 

F"  T"  V 


d'  d 


k 

2H  T" 


^  d  / 

x(y-x) 


^  ^  (y-x)' 


"Jn  t  t 


d' 


F"  F 


2 

iJL 


(y-x)’  y  A 


6e. 


x(y-x) 


4 

I_ 

I  c 

12e  ® 


3 

ZJL 


2e. 


I®H<  l«yl 


6e 


I  2 


-  M  (y-x)  -  0 


end  ®  (5  ^y 


L  (A2.10) 
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The  initial  values  (x^,  j^)  for  starting  the  iterations 
of  equations  (a2.5)  or  (A2.6)  are  obtained  as  described  in 
Appendix  I  (i.e.,  x^  =  eei  and  =  ee^),  while  those 
for  starting  the  iterations  of  equations  (A2.6),  (A2.7) 
or  (A2.9),  (A2,10),  (A2.11)  are  obtained  from  the  solutions 
ol'  equations  (A2.5)  or  (A2.8),  respectively.  If  M  is 
negative,  the  terras  e*  and  ei  are  Interchanged  (I.e., 
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X  =  and  y  =  ei)  and  the  same  equations  then  are 
employed  along  with  the  absolute  value  of  M. 

A2.3  MODIFICATION  OF  THE  METHOD 

The  method  described  above  has  been  found  to  be  adequate 
for  most  of  the  cases.  However,  It  has  to  be  slightly 
modified  in  two  cases.  Both  of  these  cases  occur  when  the 
load  on  the  arch  Is  near  to  the  ultimate  load.  As 
described  in  §2.3,5,  the  method  of  obtaining  the  ultimate 
load  consists  of  successively  increasing  the  load  by  a 
small  amount  and  checking  for  failure  at  each  step. 

The  following  two  types  of  difficulties  can  occur  during 
this  process. 

a)  If  the  load  on  the  arch  happens  to  exceed 
the  ultimate  load,  the  internal  forces  N  and  M  acting 
at  the  critical  section  (i.e.,  the  section  at  which 
the  maximum  concrete  strain  occurs)  become  too  large 
and  the  solutions  of  the  governing,  force-strain  relations 
converge  on  erroneous  values  after  going  through  a  large 
number  of  iterations.  In  a  normal  case,  it  takes  a 
maximum  of  only  fifteen  iterations  to  obtain  the  con¬ 
vergence.  Therefore,  if  for  any  case  it  takes  more  than 
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say  thirty  iterations,  to  obtain  convergence  it  is  certain 
that  the  internal  forces  are  too  large,  caused  by  too 
large  a  load.  Consequently,  the  occurence  of  thirty 
Iterations  or  more  is  taken  as  a  criterion  to  identify 
the  case  of  too  large  a  load  and  to  stop  further 
iterations . 

b)  In  some  cases  when  the  load  on  the  arch  is 

slightly  less  than  the  ultimate,  the  functions  f(x, y)  =  0 

and  g(x, y)  =  0  describing  the  force-strain  relations 

exhibit  a  peculiar  behavior.  For  clarity,  this  phenomenon 

is  discussed  here  for  a  function  having  only  one  variable, 

say,  F  CZ)=0.  The  recurrence  formula  for  the  ’Newton- 

( 12) 

Raphson  iteration’  in  this  case  is  , 


Az 


z 

k  +  i 


If  the  curve  F{  Z  )  =  2L  is  as  shown  in  Fig.  A2.1  and  we 
are  interested  in  a  solution  :z.  -  Y  ,  we  must  start 

the  iteration  at  a  value  such  that  F  (  Z  ^  )  <^  O 

/ 

If  instead,  we  happen  to  start  the  Iteration  at  y  , 

the  solution  will  converge  on  some  value  of  z  ^  y  ^ 
as  shown  in  Fig.  A2.1. 
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Figure  A2 . 1 


o  'f 

Function  F 


In  the  present  problem,  under  the  action  of 
loads  which  are  near  ultimate,  steel  at  many  sections 
yields  and  equations  (A2.6),  (A2.7),  (A2.9),  (A2.10)  or 
(A2.11)  have  to  be  used  for  such  cases.  As  explained 
before,  the  Initial  values  for  starting  the  iteration 
for  the  first  two  equations  are  obtained  from  the 
solution  of  equation  (A2.5)  while  those  for  the  last 
three  equations  are  obtained  from  the  solution  of 
equation  (A2.8).  Both  equations  (A2.5)  and  (A2.6) 
consider  no  yielding  of  steel  and  hence  the  strains 
ei  and  e^,  obtained  from  these  equations  will  be 
necessarily  smaller  than  those  obtained  from  solutions 
of  equations  (A2.6),  (A2.7)>  (A2.9),  (A2.10)  or  (A2.11). 
This  means  that  the  Initial  values  used  for  starting 
the  iteration  are  much  smaller  than  the  final  values; 
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that  is,  a  case  similar  to  the  one  when  Z  /  (F’ig.A2.l). 

The  situation  is  remedied  as  follows.  It  is 
seen  from  Pig.  A2.1  that  unless  ,  the  value  of  Z 

at  any  iteration  is  positive,  a  convergence  on  Z  =  ■/ 
is  not  possible.  If  at  the  end  of  any  iteration 
turns  out  to  be  negative  and  la.rger  in  magnitude  than  Z^, 
Instead  of  using  ^^k  [equation  (A2.12)] 

we  use. 


Z  n  z  -  Z  (AZ  J) 

k  +  l  K  k 


(A2.U) 


This  alteration  in  the  method  assures  that  Z^  will  always 
be  positive  and  a  convergence  on  Z  =  /  may  be  possible. 

The  above-mentioned  technique  is  used  to 
maintain  the  proper  signs  of  the  unknowns  x  and  y  (i.e., 
the  strains  ei  and  e^,  in  the  case  of  a  positive  M  or 
the  strains  64  and  ei  in  the  case  of  a  negative  M)  by 
using  formulae  similar  to  equation  (A2.13).  The  proper 
signs  of  the  strains  are  known  from  the  considerations 
presented  in  Appendix  I. 
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NPT1=NPT2+1 

DO  15  I=1»NSEG 
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Q0=V0*C0SF ( PH  10 ) -HO* S INF ( PHIO) 

TN00=V00*SINF( PHI0)+H00*C0SF( PHIO) 

G00=H00*SINF(PHI0)-V00*C0SF(PHI0) 

RETURN 

END 


LIST 

LABEL 

SUBROUTINE  EQEQNS 

DIMENSION  PH(31).PV(31) .PLN<  31 ).W(31)»U(31)*CK(31) .FRC(62 ) 
1FRCS(93) .EE1(31) »EE4(311 »02 ( 62 ) »E1 { 3 1 ) . 
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A5=A3*FRCS( J-l)+A4*FRCS(J-2)-2.*PLN( I ) 
A6=-A4*FRCS(  J-1)+A3^*-FRCS(U-2)+C3*PLT(  I  ) 
FRCS  (  J+1  )  =  (  A1*A6+A2'*-A5  )  /  (  Al**2+A2**2  ) 


FRCS ( J+2 ) = ( A1*A5-A2*A6 ) / ( Al**2+A2**2 ) 

FRCS( J+3 ) =FRCS( J)+C4*FRCS( j-1 )+C4*FRCS( J+2 ) 
-C5*FRCS( J+1 ) +C5*FRCS( J-2 ) 
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A6= ( A1**2+4.*A2 ) 

A3=SQRTF(A6) 

DK=(-Al+A3)/2. 

A4=( (ERTO-1. )*P*B*D*<DK-DCVR)*(D-2.*DCVR) )/(2.*DK) 
A5=B*DK*( 3.* ( D-DCVR) -DK) /6. 


SGMC=FRCA( J)/(A4+A5) 
STR1=SGMC/EBRC 
STR2=-STR1*( D-DK) /DK 
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APPENDIX  IV 


CONCISED  FLOW-CHART  AND  DIGITAL  COMPUTER  PROGRAM 
DYNAMIC  CASE 


MA I N  PROGRAM  S  UBROU  T I NES 


-214- 


DYNAMIC  LOADING  ON  ARCHES 


u 


— - 

•> 

H- 

X 

•4 

CO 

4k 

X 

_l 

2 

•> 

Cl 

Ql 

CO 

*4 

CO 

rH 

o 

Z 

•4 

» 

u 

u 

o 

» 

CO 

OJ 

_J 

H- 

Z 

o 

O 

a 

< 

3 

X 

H- 

— 

cn 

<NJ 

CO 

rH 

•4 

z 

CL 

0k 

LU 

CO 

h— 

h- 

•> 

CNJ 

a 

a. 

3 

o 

a 

*4 

01 

LU 

_i 

z 

CO 

X 

z 

CO 

1- 

#4 

h- 

— * 

CO 

a 

•k 

CO 

X 

•4* 

«4 

CL 

X 

o 

LU 

UJ 

CO 

•> 

h- 

r-i 

LU 

CO 

01 

O 

3 

CO 

•« 

o 

CO 

z 

h- 

CO 

h~ 

O 

rn 

CNJ 

— » 

CO 

z 

* 

3 

UJ 

0k 

w 

O 

CM 

o 

w 

*v 

< 

•> 

#4 

X 

h- 

UJ 

CO 

vt 

u 

< 

CM 

CL 

u 

•> 

< 

h- 

J— 

H- 

LU 

O 

Q_ 

CO 

z 

•• 

3 

X 

#4 

H~ 

CO 

U) 

. — 

LU 

f- 

vO 

r*H 

CL 

■«— 

o 

•4 

« 

o 

•> 

U 

UJ 

*4 

CO 

r-J 

CO 

u 

< 

•4 

♦ 

X 

o 

•> 

CO 

o 

W 

Q_ 

f— 

in 

X 

CL 

O 

h- 

• — 

CO 

CO 

CO 

X 

u 

1— 

O 

— 

3^ 

r-H 

< 

•> 

h- 

5 

Ql 

Ci 

CM 

UJ 

0> 

Cj 

«« 

4- 

CNJ 

CL 

-J 

z 

Q 

•> 

u 

_J 

(N 

u 

Q 

0k 

f— 

-J 

o 

>- 

r-i 

* 

#4 

CO 

4k 

a 

CO 

CO 

CO 

X 

CO 

X 

UJ 

«4 

CO 

«4 

fH 

— 

CL 

u 

1— 

UJ 

o 

#4 

4—4 

CO 

CNJ 

h— 

#> 

•« 

3 

*4 

c<^ 

o 

O 

CO 

•— 

< 

CO 

u 

CNJ 

CM 

CM 

•w 

cn 

CO 

w 

3 

LU 

iU 

u 

O 

UJ 

t/) 

w 

w 

*• 

•> 

«4 

LU 

QC 

z 

»— 

UJ 

o 

cO 

1- 

•4 

H- 

u 

3 

•> 

rH 

CO 

0. 

U 

X 

#. 

a 

<N 

CO 

— » 

«« 

a: 

* 

CM 

o 

1— 

Q 

o 

w 

o 

a 

a 

O 

«4 

c 

CO 

2 

< 

r 

cO 

> 

«4 

'd*. 

CL 

cn 

w 

rQ 

01 

LO 

< 

U 

w 

o 

CO 

is 

A 

— * 

#> 

cO 

a. 

3 

•s 

z 

f- 

v: 

CO 

a 

>- 

LU 

CM 

CO 

I-H 

3 

o 

3 

Q. 

w 

O 

X 

I—I 

0> 

«> 

cO 

«4 

O 

< 

CL 

2 

u 

3 

r-» 

r-4 

o 

z 

CO 

u 

CO 

rH 

o 

a 

Q 

#4 

4k 

3 

»4 

*4 

CO 

LO 

o 

3 

3 

1— 

Q 

4k 

CO 

r-4 

UJ 

3 

J- 

CNJ 

LU  a 

o 

z 

z 

3 

3 

3 

CL 

-J 

h- 

z 

w 

X 

» 

*2 

< 

0k 

h- 

UJ 

0k 

CO 

3 

cO 

4k 

h- 

3: 

o 

CL 

> 

a 

«4 

01 

+-4  u  CO 

CO 

4-H 

i-H 

Z 

Z 

Z 

o 

CO 

«4 

0 

c/) 

0k 

rH 

X 

a  a 

h- 

X 

II 

II 

0k 

► 

•k 

o 

— * 

•s 

> 

CO 

o 

Q 

•4 

h- 

< 

u 

o 

•4  (/) 

#4 

r-4 

1—4 

1—4 

rH 

CO 

cn 

o 

3 

4-H 

u 

QC 

h- 

o 

3 

*4 

z 

Ql 

o  - 

o 

»4 

II 

II 

II 

3 

a 

CO 

CO 

x: 

cci 

•> 

rH 

•> 

» 

u> 

a 

01 

0- 

X  01  12 

< 

— 

*-H 

1—4 

k— 4 

*— « 

3  V 

X 

w 

CO  q: 

Z 

01 

0- 

1— 

CL 

CC  O 

> 

o 

k—i 

1—4 

*4 

0k 

r— 

1— 

X  a 

z 

o 

u 

X 

CO 

w 

> 

o 

0. 

> 

X 

> 

•‘3  0 

3 

w 

w 

— 

— 

vO 

oc 

LJ  C 

CO 

2 

CO 

VJ 

_J 

UJ 

X 

U 

a  >• 

0k 

f- 

z 

X 

X 

k-H 

)-H 

k— 1 

• 

LU 

a  q: 

QC 

'W 

— 

CO 

i-H 

LU 

Q 

J— 

•< 

3 

«4 

u 

Q 

•>.  {/) 

>  4/1 

»>. 

a 

(L 

w 

w 

CM 

^  CL 

3  a 

h- 

h“ 

o 

LU 

CO  UJ 

4k 

Q 

o 

0k 

o  •> 

1— 

CO 

CL 

o 

13 

CO 

Q  3 

■-I  Cfl  O 

e? 

•« 

3 

(N| 

•s 

w 

4k 

Q 

z 

U 

O 

LU 

3 

1— 

O 

r-4  k- 

z 

X 

Ql 

h- 

CL 

3  O 

3 

*-4  CL 

z 

— 

Q 

•> 

o 

9< 

cO 

01 

01 

•4 

X 

3 

•4 

I  4/1 

r-i 

u 

«3 

(L 

CL 

h- 

1- 

1— 

4-  3 

•  o 

o 

•4 

o 

f—i 

UJ 

CO 

CO 

Q1 

Q_ 

> 

o 

X  Q. 

CO 

a  lu 

H- 

CO 

Z 

w 

w 

4— ' 

■4— 

z 

CO 

— » 

CO 

CO 

•> 

h- 

CO 

u 

U 

Z 

*4  « 

•4 

•4 

#4 

«4 

«4 

•4 

•4 

3 

•to 

w 

CO 

— > 

1—4 

z 

«> 

«4 

2 

*— 4 

rH 

1—4  rH 

rH 

4— { 

CM 

rH 

rH 

rH 

1—4 

r-4 

t- 

1-  I- 

>-  P' 

» 

X 

X 

o 

CO 

o 

X 

CO 

h- 

OD 

Q 

<  <  3 

<  UJ 

a 

»- 

2 

CO 

w 

2 

CL 

U 

z 

0k 

< 

o  o 

o  o 

Q 

o 

Q 

13 

o 

o 

2  a 

3  3; 

CO 

2 

(N 

w 

CM  2 

h- 

Ql 

Z 

o  <  < 

<  < 

< 

< 

< 

< 

< 

< 

< 

< 

(2  12  2 

-1 

w 

CM 

LU  U) 

UJ 

o 

3 

a 

h- 

^  UJ  LU 

UJ  LU 

UJ 

UJ 

UJ 

UJ 

LU 

3 

3 

3 

O  O  o 

>  o 

a 

2 

o 

•4 

U  UJ 

u 

a 

UJ 

'«k 

Z  01 

CL 

01  01 

01 

01 

□C 

01 

CL 

01 

01 

CL 

3  3  U 

rH 

CM 

co 

•+ 

»n  vO 

1—4 

CNI 

CO 

m 

w  ^  ^  ^  w  m  {O 


t/? 

z 


KD  I 


♦ 


O 

Z 

u  u  u 


in 


iH  fvi 


u 


-215- 


EPSYS=SYLDS/EST 

SPRCS=0.85*SPC 

EC= 1800, +460. *SPRCS 

EPRCS=2.*SPRCS/EC 

PHI2=PHI0+PHI0 


THTA  =  P-HI2  /  (  SEGN*PHIO) 

THTAA=PHI2/SEGN 

NSEG3=3*NSEG 


(N 

* 

O 

X 

a 

I 

og 

* 


(M 


fM 
1 

a:  V. 
X  o  ' 
u  q: 

z:  II 


X 

>!< 

O 

>K 

« 

lH 

og 

4* 


>• 

oO 

a 

* 


(/)  h- 

u  < 

CL  t- 
LU  ♦ 
O 
U  > 
I-  iA 
D  Q. 
UJ  LJ 
(L  ^  ^ 
Q 

O  O  O 
»-«  _J  ^ 
X  >-  >- 


9k 

!l 

i-H 

•k 

to 

9k 

11 

eg 

M 

• 

-J 

•k 

H* 

vO 

>- 

X 

to 

Li_ 

to 

CO 

to 

•k 

X 

«k 

X 

1—4 

n 

vO 

vO 

4 

«— 

o 

9k 

• 

«k 

u 

(X 

X 

o 

X 

Ci 

X 

CO 

u_ 

CO 

h- 

CO 

9k 

•k 

•k 

9k 

•k 

vO 

II 

o 

X 

• 

Q 

9 

CO 

CO 

Q 

CO 

•k 

Li_ 

-J 

IL. 

t—t 

CO 

9k 

>- 

■k 

— 

• 

ii 

to 

eg 

II 

o 

lA 

u 

X 

1—4 

u 

-J 

U- 

Q_ 

vO 

* 

to 

1— 

9k 

to 

«k 

II 

X 

9k 

H 

X 

X 

H 

CO 

►—4 

X 

•4* 

(O 

-J 

9k 

•k 

> 

•k 

X 

X 

U 

X 

gD 

X 

(O 

t— 4 

O 

CO 

• 

4 

•k 

X 

'  9k 

CO 

9k 

to 

uA 

vO 

Ii. 

X 

w 

cc: 

• 

vO 

CO 

CO 

q: 

>- 

X 

eg 

• 

11 

•k 

U- 

o 

CO 

i-H 

CO 

a 

CNJ 

9k 

2 

1—4 

Lu 

u_ 

u 

•— t  ^ 

.  II 

9k 

eg 

9k 

«k 

cc 

9k 

o 

a 

• 

Ii 

II 

CL 

II 

X 

i—4 

lA 

1— 

to 

UJ 

o 

> 

a 

w 

Li. 

to 

u 

X 

UJ 

u 

1— 

«k 

LU 

QC 

4 

to 

X 

to 

Q. 

II 

CO 

X 

CL 

9. 

• 

X 

lA 

H- 

O 

O 

• 

4* 

UJ 

X 

CO 

X 

9k 

•> 

1— 

x 

kO 

9k 

x 

CO 

U- 

lO 

X 

1— 

a. 

eg 

u. 

X 

vO 

9k 

9‘ 

CO 

to 

9- 

9, 

9> 

CO 

9k 

o 

II 

X 

*k 

to 

»-4 

X 

l{ 

9k 

X 

• 

h- 

CO 

eg 

X 

to 

to 

1— 

#• 

»k 

CO 

X 

vO 

CO 

CO 

< 

9‘ 

* 

AJ 

«k 

9k 

-J 

9k 

X 

• 

4^ 

X 

X 

eg 

>!< 

* 

o 

u 

o 

to 

— 1 

eg 

eg 

x> 

gD 

4^ 

i— 

H-4 

O' 

=!< 

cc 

X 

tj 

#k 

9k 

• 

9k 

X 

• 

II 

X 

*k 

1— 4 

o 

9k 

to 

X 

o 

o 

>4^ 

lA 

X 

9- 

CO 

Q 

CA 

II 

X 

cc 

cc 

•k 

X 

UJ 

X 

X 

X 

CO 

II 

X 

u 

rk 

X 

Ql 

> 

1- 

UJ 

to 

> 

o 

9k 

9k 

o 

rk 

X 

X 

X 

* 

?!< 

to 

to 

1— 

to 

♦k 

X 

X 

II 

o 

X 

II 

X 

CO 

u 

lA 

to. 

to 

C 

>- 

< 

o 

UJ 

Q 

•k 

9k 

X 

x) 

9 

X 

u 

to 

»k 

X 

• 

eg 

u 

u 

_) 

to 

X 

►- 

•k 

*k 

u 

X 

H- 

9k 

X 

CO 

X 

f— 

X 

vO 

X 

rH 

ft 

cc 

(X 

>- 

a 

h- 

X 

Q 

o 

X 

X 

X 

1—4 

eg 

X 

4 

3 

vO 

• 

lA 

+ 

eg 

* 

CL 

a 

to 

LU 

•v. 

X 

#k 

X 

X 

u 

4-4 

9k 

9k 

9k 

X 

•k 

CO 

•k 

CM 

Vs 

o 

to 

LU 

* 

* 

> 

9k 

CD 

cc 

to 

X 

1- 

X 

II 

X 

X 

X 

X 

X 

* 

a 

*— « 

* 

* 

h- 

1— 

h- 

l—i 

9k 

9k 

#k 

rk 

•k 

CO 

o 

CO 

kA 

CO 

•k 

CO 

* 

'W 

± 

U 

o 

to 

to 

X 

II 

O 

o 

LU 

o 

eg 

CO 

4 

lA 

w 

1—4 

— - 

9k 

■w 

il 

»— < 

u_ 

a 

to 

to 

to 

to 

♦—1 

1— 

►-H 

o 

o 

ii 

II 

O 

CM 

CM 

CM 

eg 

CM 

eg 

X 

X 

Q 

a 

1— 

— > 

II 

II 

II 

II 

h“ 

_J 

II 

o 

II 

X 

}— 

X 

1— 

CO 

h- 

>- 

1- 

CV 

II 

Q 

Q 

Q 

CJ 

II 

?■ 

II 

«“< 

f-H 

1—4 

I— 

1— 

1— 

f- 

1— 

1— 

< 

X 

< 

•« 

< 

to 

< 

II 

O 

h- 

U 

Q 

>- 

h- 

Ii 

»— « 

4—4 

9k 

•k 

1- 

X 

X 

X 

X 

X 

X 

e? 

X 

lA 

51 

o 

s: 

X 

x 

> 

to 

< 

IX 

X 

-J 

to 

X 

»— 1 

»— 1 

4—4 

X 

4—4 

4—4 

4-H 

•—4 

)— 4 

4—4 

LU 

X 

9k 

X 

9 

X 

X 

X 

to 

Ii 

X 

a 

X 

>- 

X 

1— t 

X 

o 

X 

2: 

w 

o 

X 

X 

X 

X 

X 

X 

to 

o 

X 

o 

ON 

o 

X 

o 

0. 

X 

1— 

to 

LU 

to 

LU 

h- 

s 

o 

h- 

o 

CD 

o 

3 

u 

X 

X 

X 

r];_ 

X 

X 

X 

X 

CO 

X 

X 

X 

vO 

X 

o 

og 


og 

rg 


eg  eg 


ir> 

eg 


-216~ 


1E8.6,3X.5HTSTP=.F8.4,3X»6HPINCR=.F6.3 ) 

FORMAT  ( 3X»6HPRDPK(  » I  2 . 2ri )  =  , F9 .6  » 3X . 6HPTGPK ( , 12 .2H )  =  . F9 . 6 . 3X . 4HTRS 
1(  »I2*2H)  =  ,F9.6,3X»4HTLD( , I  2 » 2H )  =  , F9. 6 . 3X »4HTDC (  ,  I  2 » 2H )  =  , F9 .6 ) 


NONDIMENSIONALIZATION  OF  LOADS 
DO  30  I=«1,NSEG 
PLNPKt I )=PRDPK( I )/(SPRCS*D) 
PLTPK(I)=PTGPK(I ) /(SPRCS*D) 
TRSN( I )=TRS{ I ) /TNAT 


KJ 


O 

Ctl 

Ll_ 


to 

Q 

- 

Z 

o 

u 

LU 

to 

ovi 

* 

-J 

♦ 

_J 

H- 

hh 

< 

Z 

X 

1— 

<}• 

* 

C\J 

•» 

• 

X 

rg 

< 

< 

• 

00 

X 

Lj_ 

a 

♦ 


<1- 

•  < 

(N  < 
V.  h- 

<  X 

<  H- 

h-  W 

■  J-.  Li_ 
W) 


tx 

LU 


* 

Cl 

* 

U 

Cs^ 

Cl 

to 


< 

< 

(N 

o 

o 

z 

to 

X 

Li_ 

VJ 

"S, 

h« 

h- 

1- 

< 

Z 

* 

<  X 

X 

X 

z 

< 

•—1 

< 

<  a 

1— 

< 

1— ■ 

to 

1“ 

K  * 

< 

»— 1 

X 

X 

X 

X  o 

1— 

Z 

w 

— • 

to 

I- 

< 

h- 

l~  X 

a 

1— 

o 

u 

z 

1— 

—  * 

z 

o 

o 

X 

_J 

Q 

o 

u. 

X 

o 

X  o 

• 

• 

a 

f- 

h- 

u 

1- 

1— 

t— 1 

Z  X 

o 

o 

h- 

li 

II 

LU 

O 

* 

X 

<  X 

11 

II 

II 

LU  to 

X 

LU  o 

u 

o 

a 

1-  * 

t— * 

3  I- 

•-H 

t-H 

z 

f- 

X 

* 

♦  a 

»— • 

»— 1 

Z  II 

<w> 

D 

X 

1^ 

X 

o 

* 

o 

w 

w 

•-  Z 

z 

z 

h- 

CL 

CD 

u 

X 

X 

U  CD 

X 

X 

1-  a 

Q 

u 

z 

X 

II 

il 

II 

II 

II  fi 

r 

T. 

Z  L- 

_l 

Q 

o 

o 

»— ♦ 

C\J 

cn 

m  vo 

O 

CM 

cn 

O  to 

1- 

u 

KJ 

u 

U  tJ  u 

O 

LU 

LU 

tu  1- 

c\j 

o 


lA 

vO 


in 


O 

Z 


Q 

X 

U 

(X 

< 

X 


o 


u 


o 


o 

in 


o  m 

li  m 


O 
.x  o 


o 

O 

z 

L~  CNJ 

il 

II 

Q 

X 

a 

tO 

o 

O 

•-* 

Z 

1- 

h- 

1—  X 

11 

II 

X 

Q 

o 

z 

to 

•>  CO 

1— 

1— 

< 

3 

h- 

(O 

5: 

3 

3 

O 

O 

X 

1 

0 

1— 

•> 

3 

3 

a 

2: 

1- 

3 

3 

»-4 

1— 

H-  < 

w 

w 

3 

3 

3 

Z  X 

h- 

I— 

3 

3 

3 

HH  X 

3 

< 

< 

< 

3X0 

3 

jS 

ID  O 

U 

U 

»— 4 

a  3 

m 

0 

CNJ  CO 

m 

vO 

vO  'O 

e? 

LU 

to 


•  il 

h-  HH 

IX 

<  <1- 
I—  sO 
to 

o 


vO 


217. 


PRDPK( I )=PIMCR*PRDPK( I ) 
PTGPK( I )=PINCR*PTGPK( I ) 
CONTINUE 


fH  fH 

1  1 

QC 

K 

u 

0>  * 

• 

CM 

w  w 

w 

X  X 

V, 

h~  h- 

^  vO 

3  3: 

—  U 

CM  CM 

5>" 

o  o 

h-  ^ 

n 

<-«v 

>: 

rH  h- 

rH  rH 

1  “) 

1  ( 

l>-4  «» 

X  0^ 

*— »  fH 

1—  t—  CM  CM 

h-  1 

^  X  00  CO 

_J  -H 

•k  •«  *  * 

a. 

*— 1  t— ■  <f-  m 

+  :z 

w  ^  CO  00 

^  -j 

X  X 

:i:  a 

K-  m 

H-  + 

3  3  00  CO 

~) 

CM  <M 

o 

Q  Q - 

a: 

o 

1  1  U) 

1- 

O  r-l 

-r  “) 

^  ^  a:  cn 

o 

H*  •• 

q:  rt:  >  > 

»> 

rH  O 

_J  HH 

h-  h-  u  u 

r 

a> 

CL 

^  1  I 

oc 

h- 


o 

iH 

w 

Z 

•> 

.—H 

X 

X 

3 

0 

0 

rH 

r- 

3 

»— < 

H- 

1— 

t— 

• 

• 

iH 

«« 

in 

CO 

Q_ 

w 

1— » 

r— » 

c? 

3 

Mi 

0 

0 

0 

0 

•k 

o 

i"^ 

L9 

u 

w 

0 

X 

X 

— 

3 

— * 

CM 

CM 

—• 

r 

• 

II 

It 

CM 

o 

LU 

II 

II 

3 

1— 

t~ 

3 

3 

CO 

3 

3 

0 

3 

1— 

0 

0 

rH 

fH 

CO 

00 

3 

Mi 

Mi 

K 

11 

II 

1— 

3 

II 

It 

X 

X 

iH 

CM 

z  cc: 

Q 

Z 

CM 

CM 

X 

X 

0 

Z 

— * 

h- 

1— 

3 

h- 

1 

«k 

1- 

1— 

^-l 

Z 

•s 

Q 

Q 

u 

u 

rH 

3 

Z  CM 

Vr- 

X 

X 

X 

3 

3 

r-. 

1— 

2 

i-H 

«H' 

CM 

Q 

Q 

CM 

a 

fH 

0 

CM 

m 

3 

Z 

«3 

0 

0 

11 

H-» 

1— 

1— 

fH 

X 

h- 

CM 

Z 

+ 

+ 

1 

< 

< 

II 

* 

UJ 

CO 

+ 

3 

II 

H 

II 

3 

3 

3 

0 

CO 

3 

3 

< 

3 

*— * 

3 

3 

0 

3 

3 

h- 

»— 

1 

h- 

o 

0 

•—< 

t-H 

f-H 

3 

•-H 

cc 

0 

*-4 

</) 

iO 

3 

iH 

*— r 

1— 

0 

0 

0 

X 

•> 

»> 

X 

II 

*k 

•k 

3 

3 

X 

3 

+ 

h- 

1— 

_l 

_l 

w 

Q 

h- 

3 

»— « 

*-H 

CD 

CD 

Z 

0 

CO 

3 

3 

3 

0 

0 

•— • 

1-^ 

0 

rH 

rH 

Z 

Z 

1— 

3 

?* 

3 

s: 

3 

2: 

0 

X 

X 

v: 

m 

w 

< 

< 

>— « 

0 

iH 

II 

II 

h- 

w 

3 

»- 

II 

II 

M 

_J 

CO 

h- 

1— 

h- 

3 

It 

3 

00 

3 

3: 

w 

1— 

h- 

3 

3 

3 

3 

3 

rH 

h- 

H- 

3 

1— 

1— 

h> 

■*— ' 

3 

II 

QC 

3 

■z. 

3 

cr 

3 

Mi 

Mi 

2: 

3 

3 

3 

3 

3 

*— • 

3 

yz 

rH 

Z 

3 

?: 

3 

-*£ 

3 

? 

i— 

1— 

h- 

U- 

U- 

< 

0  CM 

CM 

0 

< 

t-  < 

0 

X 

X 

3 

3 

0 

0  < 

< 

< 

< 

< 

0 

' — ■ 

CM 

CM 

CM 

CM 

CM 

CM 

3 

< 

H- 

3 

i/ 

r~r 

U 

Q 

Q 

o- 

■u 

0 

u 

U 

C3  0 

U 

d 

0 

Q 

Q 

0  0 

a 

Q 

Q 

0 

U 

in  o  LH  o  CM 

vO  00  CO 


o 

O 

00  00  r-l 


iTl 

o 


-218- 


GO  TO  115 

112  IF  (MTM-MPI)  115.114.115 


114  CALL  PRINTR 
MPI=MPI+MLT 

115  DO  140  1=1. NPT 

IF  (EKI.JTM))  117.117.120 
117  IF  (E4(I.JTM))  140.140.125 
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DWT( I .JTM-2)=DWT( I.JTM-1) 
DWT( I.JTM-1)=DWT( I.JTM) 

DWT( I .JTM)=0.0 

D2UT( I . JTM-2)=D2UT( I .JTM-1) 

D2UT( I .JTM-1 )=D2UT( I.JTM) 
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3  72  PLN(  I  .  JTM')  =PLNPK(  I  )*(TDC.N(  I  )-TM)/(TDCN(  I  )-TLDN(  I  )  ) 
PLT( I .JTM)=PLTPK(  I )*(TDCN( I )-TM)/(TDCN(  I )-TLDN(  I  )  ) 


GO  TO  400 
380  PLN(I»JTM) 
PLT ( I . JTM ) 
400  CONTINUE 
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